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Abstract 


This  report  describes  a stage  of  development  of  Surface  Ray 
Methods  for  mutual  coupling  in  conformal  arrays  of  aperture  elements  on 
conducting  cylindrical  surfaces,  with  a subsequent  generalization  to 
conical  surfaces.  Asymptotic  techniques  are  used  to  develop  expressions 
for  the  dyadic  Green's  function  for  a circumferential  and  axial  magnetic 
current  sources.  These  expressions  are  to  lowest  order  in  and 

to  two  asymptotic  orders  in  , where  ka  is  the  cylinder  circumference 

in  free  space  wavelengths  and  kD  is  the  length  in  radians  of  the  geodesic 
distance  from  the  source  to  the  observation  point  on  the  cylinder  surface. 

The  region  of  validity  of  these  expressions  does  not  include  the  paraxial 
region.  Rigorous  transition  functions  for  this  region  have  not  yet  been 
derived,  but  instead,  in  this  report,  approximate  transition  functions 
have  been  developed.  These  transition  functions  are  employed  to  obtain 
numerical  results  for  the  mutual  admittance-or  the  coupling  coefficients  between 
"single  mode"  circumferential  and  axial  slots  on  cylinders  with  kasslO  or  20. 
The  numerical  results  are  compared  with  those  obtained  via  harmonic  series 
and  via  GTD  formulation. 
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I.  Introduction 


The  theoretical  determination  of  mutual  admittance  coefficients 
in  arrays  of  aperture  elements  in  a conducting  surface  requires  the 
knowledge  of  the  dyadic  Green's  function  of  the  surface.  This  Green's 
function  (G.  F.  ) represents  the  magnetic  field  excited  on  the  unperforated 
surface  by  a magnetic  point -current  element  located  on  the  surface.  For 
the  case  of  separable  geometry  G.  F.  is  obtained  by  modal  techniques. 
However,  even  for  a separable  geometry  such  as  a conical  surface,  the 
modal  series  are  very  slowly  convergent.  Furthermore,  modal  analysis 
is  inapplicable  to  non- separable  geometries.  A more  efficient  technique, 
also  suitable  to  nonseparable  geometry,  of  constructing  the  dyadic  G.  F.  is 
that  of  the  Geometric  Theory  of  Diffraction  (GTD)  and  its  extensions  which 
we  may  term  as  Surface  Ray  Methods.  These  are  based,similarly  to  GTD, 
on  high  frequency  asymptotic  solutions  of  problems  in  canonical  geometries, 
such  as  a circular  cylinder  or  a sphere.  These  solutions  are  then  general- 
ized on  the  basis  of  locality  of  the  high  frequency  propagation  to  other 
geometries.  However  in  contrast  to  GTD  which  considers  only  lowest  order 
asymptotic  terms.  Surface  Ray  Methods  include  also  higher  order  terms. 

This  report  describes  a stage  of  development  of  Surface  Ray 
Methods  for  analysis  of  mutual  coupling  in  arrays  of  aperture  elements  on 
circular-cylindrical  surfaces  with  application  to  conical  surfaces. 

A comprehensive  survey  of  the  state  of  the  art  of  GTD  as  applied  to 
both  the  far  and  the  surface  fields  radiated  by  elemental  magnetic  currents 
(short  slots)  located  on  smooth  curved  conducting  surfaces  of  large  radii  of 
curvature  is  found  in  Ref.  LU  . Recently,  this  procedure  has  been 


successfully  applied,  after  certain  modifications  that  preserve 
continuity  of  polarization,  to  the  determination  of  far  field  patterns 
of  slots  on  conducting  conical  surfaces  [2]  . 

These  calculated  radiation  patterns  are  based  on  the  lowest  order 
asymptotic  results  [lj  and  furthermore  do  not  include  the  cone  tip 
scattering.  Nevertheless,  they  compare  favorably  with  the  numerical 
solutions  based  on  harmonic  series  backed  by  experimental  data 
[3]  . While  the  lowest  order  asymptotic  far  field  expressions  seem  to 
yield  sufficient  accuracy  even  for  curvatures  ka  sa  9.  5,  a different 
situation  arises  in  calculation  of  surface  fields,  which  are  required  for  the 
determination  of  mutual  admittance  coefficients. 

It  is  found  that  for  a current  element  on  a conducting  circular 
cylinder,  the  validity  of  the  lowest  order  GTD  result  is  restricted  to 
observation  points  not  too  close  to  the  source  and  in  addition  these  expres- 
sions do  not  hold  in  the  paraxial  and  the  circumferential  regions. 


In  order  to  determine  the  near  zone  fields  of  the  source,  v/hich  are 
essential  for  calculation  of  the  nearest  neighbor  mutuals,  asymptotic 
expansions  for  the  Green's  function  are  required  to  higher  order  in  , 


where  D is  the  geodesic  distance  between  the  source  and  the  observa- 
tion point. 

In  this  report  a rigorous  asymptotic  expansion  is  obtained  for  both 
circumferential  and  axial  magnetic  current  element  on  a conducting  circular 
cylinder.  The  expressions  are  to  lowest  order  in  ^ and  to  second  order 
in  . These  expressions  correctly  reproduce  the  H-plane  behavior  of 

the  Green's  functions  for  a circumferential  current  element. 
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However  these  asymptotic  expressions  based  on  Fock 
type  integrals  [4]  are  invalid  in  the  paraxial  region  of  a source.  This 
is  due  to  the  fact  that  the  Airy  function  approximation  of  the  Hankel 
function  ( z)  breaks  down  for  small  values  of  z.  As  a result  e.g. 

the  lowest  order  TM  (soft)  and  the  next  order  TE  (hard)  H.„  con- 
tributions for  a circumferential  magnetic  current  element  become  un- 
bounded as  the  observation  point  tends  to  cp  = 0,  i.  e.  to  the  cylinder 
generator  through  the  source. 

As  anticipated,  the  angular  region  of  validity  of  the  asymptotic 
expansions  increases  with  larger  values  of  ka. 

The  asymptotic  evaluation  of  surface  fields  in  paraxial  region 
(E-plane)  of  a circumferential  magnetic  current  element  will  necessitate 
a different  representation  than  that  used  in  this  report  and  merits  a 
separate  study.  Here  we  have  chosen  to  concern  ourselves  with  the  de- 
velopment of  an  approximate  "Modified  Planar  Formula"  for  the  surface 
Green's  function  which  provides  a transition  to  the  paraxial  region.  The 
Modified  Planar  Formula  which  is  cast  in  an  invariant  vector  form,  includes 
a provision  for  arbitrary  orientation  of  the  tangential  magnetic  sources. 

It  reduces  to  the  planar  dyadic  Green's  Function  to  of — i — — j,  when  the  ob- 

(kD  r 

servation  point  tends  to  the  source.  When  the  Modified  Planar  Formula  is 

specialized  to  the  case  of  H.  for  a circumferential  magnetic  current 

'■O 

element, it  will  be  referred  to  as  the  "Full  Formula". 

The  "Full  Formula"  matches  in  the  H-plane  the  exact  asymptotic 
expansion  to  second  order  in  ■ , is  finite  for  :p  = 0 and  reduces  to  the 
planar  result  near  the  source.  A variant  of  the  "Full  Formula"  is  also  con- 


sidered in  which  one  term  of  higher  order  in  has  been  neglected  and 


the  soft  polarization  term  is  discarded  for  observation  direction  close 
to  the  E-plane. 

In  contrast  to  the  case  of  circumferential  sources,  the  asymptotic 

expansion  to  o( — L_  J for  H due  to  an  axial  magnetic  current  element  and 
(kD)  Z (2) 

based  on  the  Airy  function  approximation  of  (z) , is  formally  well 

behaved  in  every  direction.  From  the  comparison  of  the 

numerical  data  for  mutual  admittance  between  axial  slots  on  a conducting 
circular  cylinder  of  ka=s9.  5 with  the  harmonic  series  results  [5],  one 
may  conclude  that  this  approximation  is  fairly  accurate  for  calcu- 
lation of  mutual  admittances  between  axial  slots  and  by  GTD  generalization, 
for  radial  slots  on  conical  surfaces. 

The  structure  of  this  report  is  as  follows:  / 

Chapter  II  develops  the  asymptotics  to  second  order  in  for  the 

surface  fields  excited  by  magnetic  point  sources  located  on  a conducting 
circular  cylinder.  In  Section  1,  the  scalar  TE  and  TM  potentials  are  con- 
structed for  the  case  of  an  infinitely  extended  conducting  circular  cylinder. 

In  Section  2,  an  integral  representation  suitable  for  development  of  the 
surface  ray  formalism  is  presented.  Section  3 is  devoted  to  the  evaluation 
of  curvature  terms.  Section  4 addresses  the  asymptotic  evaluation  of  in- 
tegrals developed  in  Section  2.  Section  5 discusses  the  approximate  formulae 
used  in  this  report  for  numerical  evoluation  of  the  mutual  admittance 
between  rectangular  slots.  Section  6 presents  a list  of  the  various  formulae. 
Chapter  III  is  devoted  to  numerical  results. 

Chapter  IV  deals  with  the  GTD  generalization  from  the  cylinder  to  a 

cone. 

Conclusions  and  discussion  are  found  in  Chapter  V. 

Appendix  A reviews  the  theory  of  the  alternative  representations  for 
scalar  potentials. 
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Appendix  B considers  the  details  of  analysis  of  the  curvature 
terms  (Fock  type  integrals). 

Appendix  C summarizes  the  derivation  of  asymptotic  series  for 
Laplace  type  integrals. 

Appendix  D deals  with  the  vector  form  of  the  planar  Green's 
function. 

Appendices  E and  F contain  computer  program  listings. 


II.  Surface  Fields  Excited  by  a Magnetic  Current  Element 
(Short  Slot)  on  a Large  Conducting  Circular  Cylinder 

For  evaluation  of  mutual  coupling  in  conformal  arrays  one  is 
particularly  interested  in  surface  magnetic  fields  excited  by  elemental 
magnetic  currents  placed  tangentially  on  the  surface  of  a perfectly  con- 
ducting, large,  circular  cylinder,  with  the  geometry  shown  in  Figure  II- 1. 


X 


Figure  II - 1 . Geometry  of  a Circular  Cylinder. 

Following  Ref.  [6,p.l97]  one  finds  for  the  total,  time  harmonic 
(e^Xt  dependent)  magnetic  field  excited  by  a magnetic  point  current 
source,  the  following  result: 


H(r,  r') 


-juueo  (Vxz)  (^'xzJS^r,  r') 


M + 


('TxVxz)  (v'x  7*x  z )!>“  (r,  r')-M, 

(1) 


where:  M denotes  the  magnetic  point  source  current  density,  e and 
\j.q  the  free  space  permittivity  and  permeability,  z the  axial  unit  vector, 
and  v'  and  r'  refer  to  the  source  coordinates. 


1.  Scalar  Potentials  tor  the  Circular  Cylinder  Problem 

The  aim  of  this  study  is  to  establish  the  high  frequency  surface 
ray  formalism  suitable  for  mutual  coupling  applications  in  cylindrical 
and  conical  arrays.  This  is  facilitated  by  the  particular  integral  repre- 
sentation of  the  potentials  S*  and  S . From  Appendix  A (9),  (13),  (19)  and 
(21a)  one  finds 

! • \Z~Z  ^ -cos/T^n-lc^-p'l) 

<**»  ? dXv — 


4TT* 


“"AT71 


r (2) 


(1) 

QH  a—  (k  a) 


7*7 


(2) 


(2) 


— lH/r-<V<> — rzj Hsr  (kt°<)JHvr;(k t0>]- 

V V AU  /l.  V V 


QH 


X ^kta) 
v 


(2) 


where  for  s'  , Q=l,  for  S*  , Q=  , and  k^  = J k2  - X , 

With  Xv  = v2  , the  contour  C,,,  in  the  Xv  -plane  maps  into  a straight 
line  from  -■  to  +®  just  below  the  real  axis  in  the  v -plane. 


s(r,r')  = 


- j nr  | z-  z'  | ® - js 

1 £ dX  e r dv  ^cosvdT-h  -gll 

4tt  “ " c W 2j  /X  J .. 

z w u -■  -)6 


sin  v rr 


V Lhv  <k,  »<  > - -nr  <V<  'JHv  <y  >>-r  • 


(1) 

QH,  (k  a)  (2) 


(2) 


QHV  (kta) 


(3) 


where  6 is  an  arbitrarily  small  positive  quantity.  The  only  singularities 
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In  order  to  cast  S into  a form  appropriate  to  the  surface  ray 
formalism,  it  is  necessary  to  convert  the  angular  Green's  function  into 
a traveling  wave  series. 


where  cpj=  |cp-cp'|,  = 2rr  - | tp -cp ' | and  -rr  < < tt  . Consequently 


2 • ®+j6  ' ®-j6 

S(£>£,)  = ' 15TT  I I J dkz  j 

i=l  4=0  -®-j 6 ' -»-j6 


-jk  Iz-z^  -jv(CP.  + 2n  4) 
dv  e z e 1 x 


(1) 

QH  (k  a)  (2) 


r wn  ye.)  -|  , 

lHvV<> m-5-  Hv  ,ktc<’j  Hv  <V>>7Z 


QHV  (kta) 


-iv®. 

One  notes  that  for  a given  value  of  v,  the  two  waves  e and 

-jvc P2 

e propagate  in  opposite  directions.  The  terms  with  4^0  represent 

waves  that  arrive  at  Cp  after  having  traveled  4 complete  circuits  around 
the  cylinder.  If  the  cylinder  radius  (ka)  is  not  too  small,  the  terms  with 
4^0  are  negligible,  because  of  the  exponential  attenuation  of  the  creeping 
rays,  as  will  be  shown.  Also  for  the  same  reason,  if  > >cp^  , the  con- 
tributions of  the  4=0,  co  ^ term  maybe  neglected  compared  to  that  of  cp^  term. 
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Thus,  it  will  be  sufficient,  except  near  = rr  , to  consider  only 


the  4 = 0 , "Pj  term,  i.  e. 


- j p -jkzU-z'|  * 

S(r,£  ) = - I dkz  e ! dv  e J 


p (1)  QH  (k  a)  (2)  (2) 

LHv  v<'— pr-Hv  <V<’JHv  <V>> 


QHV  (kta) 


In  the  above,  cpj  was  replaced  by  $ and  the  integration  contours  have  been 


shifted  to  the  real  v and  k axes. 

z 


I 


2.  Integral  Representation  for  the  Field  Components 

There  are  two  independent  tangential  source  configurations 
on  the  surface  of  a conducting  cylinder  giving  rise  to 
different  solutions.  One,  an  axial,  the  other  a circumferential 
(transverse)  magnetic  current  element. 

a.  Circumferential  Magnetic  Current  Element  (Short  Slot) 

In  this  case 


M = M cd 

— r 


and  therefore,  from  (1)  one  has 


M 


Hl  (£•£')  = s"(£*£)  ‘ ^ eo  Mcp  s'(£>£'> 


(9) 


M , ,2 

HZ  <£•£>  - I&rr  (-vt > s"  ^ • 


(10) 


where  the  superscript  "c"  denotes  quantities  excited  by  a circumferential 
magnetic  current  source. 

In  view  of  the  relation  [Ref.  6,  Eq.  33b,  p.198] 


Vt2  S"(r,r,)=.7t2s"(r,r/)  = g"  (r.r') 


/2  c" 


(ID 


(10)  may  also  be  written  in  the  form 


M .2 

D o 


TTC  / /\  P O — " / /* 

Hz  ^£*£  ) = juu w,  c'  Sz'ocp'  G ^£’£  ^ 


(12) 
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l 

J 
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The  formal  modal  representation  of  the  scalar  Green's 
function  G*  is  [Ref.  6,  Eq.  33b,  p.  198] 


G"(r,r')  =1  f"  (o)  *[V)  gz  (z.z)  . 


Upon  combining  (A2),  (A18)  and  (13)  one  finds  that  G"  satisfies 


(V2  + k2)  G"(r,  r')  = - 6 (r-r')  , 


with 

a 


G"(r,r')  | . = 0 


o = o =a 


There  is  no  TM  contribution  to  the  z-component  of  the  magnetic 


field. 


Upon  introducing  the  Wronskian 


(1)'  (2)  (1)  (2)' 

R (z)  H (z)  - H (z)  H (z)  = 


4j 


TTZ 


(13) 


(14) 


(15) 


(16) 


one  finds  that  for  o = o = a»  i*  e.  the  case  of  surface  magnetic  field  , 
excited  by  a circumferential  magnetic  current  element  located  on  the 
surface 


H*  = H"c  + H'C 

-P  =P  =P 


where  and  denote  the  TE  and  TM  contributions,  respectively. 


(17) 
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From  (8),  (9)  and  (16)  one  has 


jM  T - jk  | z-  z'\  k 

„»c  / M CP  f „ z z 

h¥  (£•£)!=  t~2 j-  J dkze  rr. 

* , 4n  uiu  a k. 

c = o -a  o t - 


2 « 

2 -iv^2Hv<kta> 

T“  J dve  v ~UT> 

t -»  H v (kfca) 


-jM  uus 

„ /C  . /.I  1 z o 

H„  ( r . r ) I = 2 — 

5 , 4TT  a 

r = o=a 


J dkz  6 


-jkjz-z'i  : 

T~  J dve 


(2)' 

-jvf  \ 


Hv  (k.a) 


Similarly,  for  , one  obtains  from  (8),  (12)  and  (16) 


H°  (r,r')  | = sgn(23-T,)sgn(z-z') 


" -jk  fz-z'l  k “ $ ^’(ka) 

■ J’-  - 2 — I dve  1 v — — 


o = c = a 


417  uu  u a 
o 


; dk  e 
J z 


Y' 


f 1,  u > 0 , 
s§n  u = 1 -1,  u < 0 . 


The  reason  for  appearance  of  the  sgn  u in  (20)  is  the  operation 

^ in  (12);  the  absence  of  sgn  u in  (18)  and  (19)  is 

^2  6 2 

brought  about  by  the  operations  an<^  ^ § ~r>  i-n  (9)  which  because 

of  the  presence  of  |z-z'|  and  Itt-a'I  in  the  exponents,  do  not  give  rise  to 
sign  changes. 
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b.  Axial  Magnetic  Current  Element  (Short  Slot) 
In  this  case 


M = M z 
— z 

and  in  view  of  (1),  one  has 


(22) 


M 7 

■ (Vx  7xz)  (-7p  S"(r,r') 


7 x V x z G*  (r,  r7)  , 


(23) 


where  the  relation  (11)  has  been  used.  Therefore, 


H*  <£•£'>  = i«nr-  ( 


M-  - s2 


3 z 


, . 2\  _ » , 

+ k j G (r  ,£ ) 


and 


M .2 

z 3 


T / / \ Z O _ H a 

H-p  <£>£  > =T¥m"T  W5T  G {£’I} 


(24) 


(25) 


Following  an  analogous  procedure  as  in  (18),  one  obtains,  via 
(8),  (11)  and  (16),  the  following  expressions 


-jM  r 

H*  (r,r')  | = — 2 — j dk  e z 

/ 4tt  ul  u a 
0 = o = a o - ® 


and 


-jk  Jz-z'l 


. (2) 

, r , -jv*  Hv  <kta) 

kt  J dVe  ~1 zr> 

Hv  (ka) 


(26) 


/ / z r 

H*  (£,_r  ) I = sgn(^-3,)*gn(*-z') — ^ Z"  ! dkze 


JM,  r ..  -jkJ2-z'lkz“  .jv*  H(Z)(k  a) 

|dve  J M 


p=  0 =a 


4H“um  a 

o -® 


z 

IT  P 


Hv''v 


(27) 


14 


where  the  superscript  "a"  denotes  the  field  excited  by  an  axial  magnetic 
current  element. 

It  is  noted,  following  a similar  reasoning  as  for  the  case  of  a cir- 
cumferential magnetic  current  element,  that  does  not  change  sign  for 
z < z or  CD  < but  may  do  so. 

3.  The  Curvature  Terms 

Relations  (18- ZO ) and  (26-27)  involve  two  contour  integrals  with 

respect  to  the  angular  and  the  axial  wavenumbers.  They  contain  two 

independent  large  parameters,  one  ka  the  circumference  of  the  cylinder  in 

wavelengths,  the  other  kD  the  free  space  phase  delay  along  the  geodesic 

distance  from  the  source  to  the  point  of  observation.  The  asymptotics  will 

be  carried  out  to  lowest  order  in  — in  each  polarization  and  to  of K~nr , 

, j \kD)3/2; 

in  the  deep  shadow,  or  equivalently  to  O1 j/in  the  near  zone.  An  evalua- 

1 *kD) 

tion  to  the  lowest  order  in  is  insufficient  for  applications  to  mutual 
coupling  in  conformal  arrays,  because  the  elements  of  the  array  are  not 
necessarily  far  apart.  Furthermore,  for  the  case  of  circumferential  slots, 
the  lowest  order  TE  contribution  vanishes  in  the  H- plane  and  a consistent 
expansion  to  second  order  in  is  essential,  as  will  be  seen  in  Section  4. 

One  notes  that,  in  principle,  it  is  sufficient  to  carry  out  the  evaluation 
of  the  fields  in  the  first  quadrant,  0 < 9 <— , those  in  other  quadrants  follow, 
in  view  of  (18-20  and  26-27)  via  the  various  sgn  u relations.  This  procedure, 
however,  will  be  replaced  by  an  equivalent  method  that  conforms  naturally  to 
the  surface  ray  description.  In  this  simple  alternative  scheme  that  avoids 
the  use  of  the  sgn  symbols  and  will  be  explained  with  the  aid  of  final  form- 
ulae, the  observation  angle  3 (Figure  II- 7 ) ranges  over  all  four  quadrants. 
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Following  a standard  procedure  (e.  g.  Ref.  4)  one  transforms 
variables  via 

v=  kta+  ftt 

where 


(28) 


(29) 


so  that 

v2  = a2  + 2kt  a ft  t + f^  t2  (30) 

Substituting  (30)  into  (18)  and  making  use  of  the  Airy  function  approxi- 
mation for  the  Hankel  functions  to  the  lowest  order  in  ^ (see  Appendix  B), 
one  finds 


tt  «ct 


_ . 4n  um  a 
P=  0 = o 


® -jk  U-z'l  - jk  af  k2  f2 

fdkze  2 4 -^VQ(x) 


“ -jk  |z-z'|  - jk  atk 

+ [dk  e Z ^ V,(x) 


kt  1 


P -jk Jz-z'l  - jk  a$  k2  fj* 

I dk  e 2 1 


3-V2(x)j 


(31) 


where  x = fj  and  V (x)  denote  the  Fock  type  integrals: 
t n 


v , , r -j*t  w2(t) 

V (x)  = I dt  e J 

0 T 

" -jxt  , w2(,) 


(x)  = i dt  e”J  t 


(32) 

(33) 


and 
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(34) 


V^(x)  = I dt  e 


-jxt  ,2 


V « 


The  integration  contour  T is  depicted  in  Figure  II - 3 . 

The  validity  of  the  Airy  function  approximation  of  Hankel  function 
H^(z)  is  predicated  upon  the  condition  of  large  z and  |v-z|  < 0 (z^^). 
If  z is  not  large  as,  e,  g.  in  the  paraxial  region  ) «0, 


Figure  II- 3 . Contour  of  Integration  for  Curvature  Terms, 
where  the  value  of  kt  a at  the  saddle  point  as  in  (31)  becomes  small,  this 
approximation  of  (kta)  breaks  down,  and  the  above  method  of  asymptotic 

evaluation  becomes  invalid.  The  extent  of  the  paraxial  region  depends  on  the 
value  of  ka;  the  larger  the  value  of  ka  the  narrower  will  be  the  paraxial 
region  where  approximation  (Bl)  does  not  hold.  To  date  no  rigorous  asymp- 
totic transition  function  is  available  for  the  paraxial  region,  but  in  Section  5 
we  shall  derive  an  approximate  transition  function,  which  will  be  employed 
in  Chapter  III  to  obtain  numerical  results  for  mutual  coupling  between  rec- 
tangular slots  on  circular  cylindrical  surfaces.  In  spite  of  this  deficiency, 
it  will  be  seen  from  comparison  of  the  numerical  results  with  those  obtained 
via  harmonic  series,  that  the  overall  agreement  is  very  reasonable  for 
ka  10  and  20,  including  the  paraxial  region  for  circumferential  and  axial 
slots. 
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where 


1 + ^ x^  + x9/2  +.  . . , for  small  x (x  < 0.  6) 

1 ^ ./T  32 


v2(x)  =, 


4v  TT 


5/2  r-  - jxt 
X . — J ] 


t e F , for  large  x (x  > 0. 6)  . (4i 


3 e 


■J  1 


4 p=l 


. TT 

-J  ^ 7 / 

Here  t = jt  | e , Ai(-lt  |)  = 0 and  Ai  (y)  denotes  the  derivative  of  the 
p p p 

Airy  function  Ai(y). 

Another  type  of  Fock  integral  that  represents  the  TM  polarization  contri- 
butions is  defined  by 


r _ixt  w2(t) 
u„(x)  = ! dt  e-Jxt 


w2(t) 


(4: 


where  t is  the  same  contour  as  in  Figure  II-3. 

Similarly  to  V (x),  U.(x)  becomes  a residue  series  for  large  values 
n u 

of  x and  may  be  expanded  in  a power  series  for  small  values  of  x as  follows: 


. TT 

4 


UQ(x)  = -2nj  -_.3T2.  Uq(x)  , 


2 y/ff  X' 


(4; 


where 


r 


1 


1 V lTT 

1 --L-* 


«„(«) 


3/2  + 4% 

for  small  x (x  < 0.  6) 


7 rz  3/Z  r - jxt 
2^x  V e P 


for  large  x (x  > 0.  6). 


(43) 


-J  T 


p=l 


. rr 
-J  T 


Here  t = !t  | e ^ ,with  Ai(- |t  |)  = 0. 

P p P 

Figures  (II- 6 , II- 8)  show  the  residue  series-power  series  crossover 


regions  for  the  curvature  terms  Vq(x),  v^(x)  and  Uq(x).  From  the  computed 
data  it  is  found  that  the  smoothest  crossover  is  obtained  for  x = 0.  6,  where 
the  difference  A between  the  sum  of  ten  residue  terms  and  the  sum  of  four 
power  series  terms  is  as  follows: 


A |vQ(x)  | = 0.  0025% 
Alv^x)  | = 0.  42% 
A|v2(x)  | = 1.7% 
a|uq(x)|  = 0.2% 


A arg  v (x)  = 0.  012° 
A arg  v^(x)  = 0.  126° 
A arg  v2(x)  = 0.  101° 
A arg  Uq (x)  = 0.  10° 


In  actual  computations,  the  crossover  value  of  x = 0.  8 was  employed,  with 


A |v  (x)  | = 0.  13% 


A arg  v.(x)  = 0.  012c 


aJv^x)  | = 0.  77%  A arg  v^x)  = 0.  13° 

a!u0(x)  | = 1.  26%  A arg  uQ(x)  = 0.  28°  . 

The  integrals  V^(x)  arise  from  the  TE  polarization  in  the  cylindrical 
geometry  or  from  the  Neumann  boundary  conditions.  In  the  terminology  bor- 
fored  from  acoustics,  the  latter  corresponds  to  a hard  boundary.  Accordingly, 
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Figure  11-4.  Curvature  Term 


I 


the  TE  polarization  is  termed  the  hard  polarization,  V^(x)  are  the  hard  curva- 
ture terms  and  the  residue  contributions  arising  from  V (x)  at  poles  t will 

n p 

give  rise  to  hard  surface  (creeping)  rays.  The  Vj(x)  and  ^(x)  are,  within  a 
numerical  factor,  the  first  and  second  derivatives  of  Vq(x).  The  TM  curvature 
term  Uq(x)  is  called  the  soft  curvature  term,  since  it  corresponds  to  the 
Dirichlet  (soft)  boundary  conditions.  The  residue  terms  resulting  from 
evaluation  of  Uq(x)  at  the  poles  t will  give  rise  to  the  soft  surface  rays. 

4.  Asymptotic  Evaluation  of  Surface  Fields. 

a.  Circumferential  Magnetic  Current  Element. 

Returning  to  (31)  and  denoting  by  D the  geodesic  distance  between 
the  source  and  observation  point  (see  Figure  II-7) 


Source 

location 


-Observation 

point 


q9>  r 


Observation 

point 

Surface  ray 


Magnitude  dipole 

Conducting 
circular  cylindar 


Figure  II-7.  A Surface  Ray  on  a Circular  Cylinder 

(a)  Developed  Cylinder,  Source  located  at  (acD',z') 

(b)  Physical  Structure  . 


D = [(z-z  ')2  + a2(cp-CD7)2]^  , 


one  finds  that  for  observation  points  in  the  first  quadrant  the  exponent  in 


(31)  may  be  rewritten  in  the  form 
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whe  re 


k \z-z'\  + k a*  = D(k  sin  9 + k.  cos  9) 
z t z t 


z-z  1 . a 

~HT$  = tan  ° ' 


Introducing  the  change  of  variable 


k = k sin  a 
z 


kfc  = k cos  a,  (4 

where  the  positive  sign  is  chosen  in  order  that  for  kz  = 0 (a  = 0),  k^  = k,  and 

\ 

substituting  (47)  and  (48)  into  (*^5),  one  finds 

k I z-z  7 1 + k.  a?  = kD  cos  (tt  - 9) , (4 

z t 

Combining  (49), and  dk  = k da, with  (31),  one  obtains, for  observation  points 

Z L 

in  the  first  quadrant , the  hard  (TE)  contributions 


H:C(r,r') 


: ■/  4tt* 

P=  p =a 


^[k2  f2  r da.'*0  «•  'a-9|co,2/5«5ln2av0(x) 

L P 

+ k2  r da,-ikD  co“  (1-S)  ,in2a  VjW 


♦ tj  ! da  e-ikD  c°a  <a-S>,inZ*  to.4,Jav2M 


where 


-(¥) 


The  integration  contour  P is  the  a-plane  is  depicted  in  Figure  II-8. 
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Top  Sheet  of  Fig.  II- 2 


Figure  II- 8 . 


One  Period  of  the  a-plane  and  Contour  of 
Integration.  Shaded  Regions  are  Valleys  for 

e-jkD  cos  a . The  Numbering  Corresponds 
to  Mapping  of  the  Four  Quadrants  of  Figure  II- 2. 


To  perform  the  asymptotic  evaluation  of  the  integrals  in  (50),  one 
first  deforms  the  contour  of  Figure  II-8  into  the  steepest  descent  path 


(SDP)  through  the  saddle  point  ag,as  depicted  in  Figure  II-9.  The  SDP 
intersects  the  real  a-axis  at  an  angle  From  (Cl)  one  has  for  the 


Im  a 


Figure  II- 9.  SDP  in  the  a-plane 

first  integral  in  (50) 

2 

q(a)  = -j  cos  (a-0)  = -j  -s 
F(a)  = sin^a  cos2^3a  Vq(x) 

and 

Of.)  = F«t>  g . 


(52) 

(53) 

(54) 
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From  (52),  the  saddle  point  is 


a = a 

s 


(55) 


The  even  power  terms  in  the  Taylor  series  expansion  of  (54)  about  s=0  are 
found  below 

G(O)  = F(0)a'  (56) 


G'(0)  = F"(0)  (a')3  + 3 F '(9) a'  a"  + F(9)  a'"  , 


s s 


(57) 


whe  re 


.(n)  _ d a 


s , n 
ds 


(58) 


s=0 


and 


F(n)(0)  _ dnF(a) 


da 


(59) 


a=9 

From  (53)  one  finds 

2 - 1 / 3a  _j_3o  tr  /..  v . -> 5/3fl_;_  n ,,  /„  \ _l  ___2/ 3_^_2( 


F '(9)  = --|cos"  ' 39  sinJ9  VQ(xs)  + 2 cos /J9sin  9 VQ(xg)  + cos  sin  9 Vg(xg)  xg  , 

(60) 


F"(9)  r --|cos'4^30  sin4  9 V0(xg)  cos'1/ 3 9 sin3 9 Vo(xs)Xg  --y- cos2/ 3 9 sin2 9 V0(xg)j 


+ 2 cos®/ 39  V (x  ) +4  cos3/30sin  9 V '(x  )x  ' + cos2^39  sin29  V^(x  )(x ')2 

US  USS  US3 


+ cos2//39  sin29  VQ(xg)  x"  , 


(61) 


where 


x = x 


kD  f 1 / 3n  X 
= f COS  0 $ 


(62) 


’a  2 i 

a=0  g 


with  f defined  in  (51)  and 
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(lie.)  ( k*  '!/3 

* = = [7Z&J 


= - x tan2  9 - i x 
9 s 3 s 


The  mapping  derivative  can  be  found  via  (52)  and  the  proper  branch 
is  determined  by  the  condition  that  at  1=0,  arg  ^ . Therefore, 


dOt  (rr  (.  .3 

HI  = ./Zj  V1  ‘ 3 T 


a # = ..^7  . 

s J 

Consequently,  the  higher  order  derivatives  become 

*:  « o. 


a » «i^3L 


As  a result,  one  has 

C(0)  = v'ST  si**2®  cos*"/3©  Vg(xg) 


G'(0)  = |cos'4/3°  sin"*  9 1 Co,_  1/3®  «in3e  Vo(x8)xs 

- ^ cos2/38  sin26  VQ(xs)  + 2 cos8/39  VQ(x9)  + 4 cos5/30  sin 9 VQ/(xg)  xg' 
+ cos2/30  sin2®  Vj(xg)(xg')2  + cos2/3  sin20  V^(xg)  xg  J 


+ sin20cos2/30  V0(xg)  . 
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Thus,  one  has,  with  0 = kD 


-jM.„k2f2| 


<!(£.£)  I ~ ~ P l ^T sin29  cos2/39  V.(x  ) 2-^2 

’ 0=  L ^rr  ° s (kD)* 


' 4*1  uuu  a 

Da  0=a  O 


+ ^v|fJJ(2c0s8/39  .^cos2/39sin2Q.|=og-4/3Qgin29)v^(x^ 

+ (4cos^^Ssin9_^cos  ^30  sin30)V  *(x  )x  ' 

3 0 s s 


2/  3 g 2q_r/,  . // 

+ cos  °sin  o V (x  )x 
O s s 


+ cos2//36  sin29  V 


o'*3X<)2J 


-jkD 


where  the  second  subscript  ' l1'  denotes  the  first  integral  in  (50). 

Carrying  out  the  various  derivatives  in  the  bracket,  one  finds 


■jM.kV 


' / 4 tt  -X  u a 

D=  0 =a  o 


VQ(xs)  cos2/30  sin29  ^ 


-jkD 


M,k2f2 


4ttSji)u  a 
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^ V0(xs)(-  |cos-4/30sin49  + 2cos8/36 


61  2/3s  . 2a N,  e "JkD 
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M k *- 

4tt2U)4  a 12 
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• v j (xg)(- 3in49  + 5cos29sin29j 


20  . 29\e'ikD 

o sin  o t 

' v£d" 


where  the  term  with  V2(Xg)  has  been  neglected  due  to  its  lower  order  in 

££■•  One  observes  that  the  three  terms  in  (73)  have  the  same  exponential 

order  in  ka,  because  the  residue  series  for  VQ(x),  V^x)  contains  the  same 

hard  exponents  (see  36  and  38).  However,  the  algebraic  orders  in  — - of  the 

ka 

terms  differ.  Therefore,  for  consistency,  the  third  term  is  also  discarded. 

If  one  transforms  VQ(x)  into  the  near  zone  form  vQ(x)  via  (35),  one 
finds  the  form 
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u'c  / I -xx  k^Y  . Zo  , , 

- - jM9-2F31a  6 v0(X3)-kD“ 

p=  c *=a 


k*TYr  2fl  61  a;2fl  lsin49vi  , v e' 

TttAcos  9 -24  sin  0.? rJv0(xg)- 


35  2tt  \ 


where  for  small  values  of  xg  one  uses  for  v0(xg)  the  power  series  expres- 
sions and  for  large  values  of  xg  the  residue  series  (see  (36)  ). 

The  second  integral  yields 


H;C,(r,r#) 


-M35k 


1 / 4rr  u)u  a 
0=  P =a  o 


UcD)‘ 


,c  . 1 /I  1/3 

Because  2 is  0(^),  it  is  also  discarded,  as  compared  to  O of 

first  two  terms  in  (73). 

The  TM  contribution  follows  via  an  analogous  procedure.  From 


(19)  one  finds 


H^r.r') 


i.^2  (x). 


1 4tt  a 

P=  P=a  a P 

where  the  integration  contour  P is  same  as  for  the  TE  case. 

The  asymptotic  evaluation  yields 

/C  / I - jkfJ  . rr=r  -jkD  " 

+ 0^ka)5''3^°(kD)3/Z  ) 

The  O'r — - Jterm  results  in  the  near  zone  in  — - — =-.  To  Cl — - — =•  J 

(k  D)3/Z  (kD)3  (k  D)2 

in  the  near  zone  one  needs  to  consider  only  the  leading  soft  contribution. 
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, 2 . -jkD 

H (r,  r ) - M uo(xs}  —7 

o=p'=a  cos  9 (kD) 


(78) 


Combining  (74)  and  (78),  one  has  the  total  x-component  of  the  surface 
magnetic  field  excited  by  a circumferential  magnetic  current  element: 


| ~ 3in20  vo(xs)  - 


-jkD 


o=p  =a 


kD 


r x.  kZY{  2 a 61  . 2fl  1 sinV]  , . 

+ M „ -= — » cos  0 - "^~r  sin  o =5— Jv  (x  ) 

qp  2*t  \ 24  9___2g/  o s 


cos 


-jkD 
0 s'  (kD) ^ 


+ M.. 


k2Y  1 


u (x  ) 


■ jkD 


cos 


(kD) 


(79) 


Besides  those  already  neglected  additional  terms  of  order 


1 


(kDp 
H<2>'(z) 


but  of  lower  order  in  ^ appear  when  the  uniform  expansion  of  ^ 

is  used.  These  terms  are  similarly  discarded. 

The  above  asymptotic  evaluation  is  valid  as  long  as  kacos  0 is  not 

too  small,  because  otherwise  the  Airy  function  approximation  breaks  down. 

(2) 

The  failure  of  the  Airy  function  approximation  of  (z)  for  small  values  of 
z manifests  itself  in  appearance  of  a spurious  algebraic  branch  point  at 
ct  ■=  ^ in  the  various  integrands  and  of  ^2and  $”^2  factors  (see  30,  35 
and  62)  with  Vq(x)  **  — - , Uq(x)  - ■ * . - —jjy  as  5-0  or  . Although 

the  proximity  of  the  saddle  point  to  the  branch  point  for  paraxial  propagation 
could  be  taken  into  account  by  appropriate  parabolic  cylinder  transition  func- 
tion  [7],  the  presence  of  and  factors  yields  infinity  in  the  axial 

direction  (§=0)  which  cannot  be  removed  by  constructing  a transition  function 
based  on  the  assumed  validity  of  (Bl).  An  approximate  transition  function  and 
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its  variant  which  render  the  field  finite  for  $=0  will  be  given  in  Section  5. 

From  (20)  one  finds,in  a similar  manner,  for  the  z-component  of  the 
surface  magnetic  field  to  lowest  order  in 


HC(r,  r ')  - sgn(co-co  ')  sgn(z-z  ') 

z , 4rr  U)|i 


3 k2‘2  f da  , -ikD  <a-  9)si»a  co.5/3a  v„(.) . 


P=P  =a 


oa  P 


(80) 


Following  the  same  procedure  as  in  (79),  one  finds 

2. 


-jkD 


H^(r,r')  j -sgnfco-a'jsgnfz-z')  j sin  9 cos  9vQ(xg) 

-jkD 


p=p  =a 


lCP  2TT  V 8 


M .^=*1  - -s-  sin  9 cos  9 + — 


5 sin^9N 


9 cos 


F J v0(xs} 


(81) 


(kD) 


Again  the  second  order  terms  are  inadequate  as  9 . As  mentioned  previ- 

ously, the  sgn  symbols  in  (81)  may  be  removed  when  the  range  of  9 is 
extended  to  all  four  quadrants  i.  e,  0 < 9 < 2^  . 
b.  Axial  Magnetic  Current  Element 

In  contrast  to  the  case  of  circumferential  magnetic  current  element, 

(2) 

in  spite  of  the  fact  that  the  Airy  function  approximation  of  H'  (z)  breaks 


(— ■ ^rJare 

\kD)2' 


down  in  the  paraxial  region,  the  asymptotic  series  terms  O' 
bounded. 

The  z-  and  cp-  components  of  the  surface  magnetic  field  may  be 
obtained  from  (26)  and  (27).  For  mutual  coupling  between  axial  slots,  one  is 
particularly  interested  in  Proceeding  in  a similar  manner  as  in  Section 

(4a)  with  (28)  and  (32)  one  has 


H 


z 


ifel-  r dae"jkDco9  ^a-S^cos8^3a  VQ(xg) 


(82) 


, 4TT  UJU  a c__. 
p=p  =a  o SDP 


where  SDP  is  shown  in  Figure  II- 9 and  VQ(xg)  is  the  hard  curvature  term, 
with  xg  given  by  (62). 
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One  finds  in  (82) 


a = 0 
s 

F(CL)  = cos8^3a  VQ(x) 

F"(a)  =^cos2/3asin2av0(x)  - 8 cos8/3a  VQ(x)  - 
+ cos8/3Vj(x)  (x')2  + cos8/3a  Vq(x)x" 


1_6 

3 


cos^3a  sintt  Vq(x)  x ' 


Therefore 

G(0)  = JIj  cos8/39  VQ(xg) 
and 


G"(0)  = 2jVir  [Trcos^e  «»26  V0(xs)  -fcos^S  VQ(*s) 

« cos5/38sin9  V '(*)*'  + cos8/30  vj(xs)(<)2  + Co. 8,3  8 


+ ii2Lcos*39V0(xs) 


8/3, 


In  accordance  with  (C6),  one  has 


H*(r,r') Z 

Z / 4 TT 
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+ cos2/39  sin2 9 -l|cos8/39)  V0(xg) 


.^cos5/39  sin  9 V(J(xg)xs/  + cos8/39  V(J(xg)x'/ 


Carrying  out  the  derivatives  in  the  bracket  yields 

-jM  kV  an  -jkD 

H*(r,r')| y JETj  V (x  )cos8/39 

Z I . A TT^illLl  a 0 S TkD 
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T^Vl(Xs)(-|COs2e  sin20+^cos40  j • -2  . (89) 


(kD) 


Note,  the  last  term  is  of  O W7T  below  the  leading  term  and  is  therefore 

(ka)2/i 

neglected. 

Finally  one  arrives  at  the  following  simpler  form 


) | VCOS  0vO(Xs) 


k2Y  2a  -jkD 

kD 


P=  P =a 


(90) 


Because  of  a similarity  between  (25)  and  (12),  component  excited 


by  an  axial  magnetic  current  element  becomes 

~ jMz^rf  sin9  cos9  v0(xs)SkF“ 


P=P  =a 


k Y/  23  . a a 5 sin  9 . . e J n^A-r-jrr 

- M -*rz-  { -~r  sin  9 cos  9 +^r n j v.(x  ) * , 0 < 9 <2n 

z 2tt  \ 4 9 cos  9/  0 s ^£^2’  — 


(91) 


5.  Approximate  Formulae 

a.  The  "Full  Formula"  for  Circumferential  Current  Element. 

Since  the  limits  of  validity  of  the  rigorous  asymptotic  expansion  (79), 
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iu 


(81)  and  (91)  do  not  extend  to  the  paraxial  region,  we  shall  derive  an 

approximate,  simple  transition  function  which  is  finite  for  9 “*  matches 

to  of — — jjthe  rigorous  formula  at  9 = 0 and  reduces  to  a planar  result  to 
(kD) 


( 1 \ 

Ol =-iwhen  the  observation  point  approaches  the  source. 

(kDj 


The  derivation  of  the  "Full  Formula"  is  given  below. 


Starting  with  (31),  we  immediately  discard  the  integral  containing 

1 


V^x),  since  it  gives  rise  to  a lower  order  of  than  considered  here. 
The  first  integral  in  (31)  is  replaced  by 

,2 


y a2  r 

2 - ? I 

4 n uuii  a o z 

o -® 


, , dk  e 
2 1 z 


-ik  Iz-z  ' | -jk  a*  f 

z x — V (x) 

kt  Vx>  ' 


while  the  second  integral  is  left  in  its  original  form  which  contains  k 

z 

Hence, 


H:C(r,r') 
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2 ' 2 
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p=  p =a  o P 


[Ji  l’dae--'kD<:os(a-0)£?V„W 


t O' 


The  contour  P is  deformed  into  SDP  of  Figure  II- 9 and  the 
asymptotic  evaluation  yields  the  leading  terms 


dae“jkDcos(a"9)f?  V„(x)  ~f2  VQ(xg)  cos2/39  /^-e_jkD 


t O' 


V kD 


SDP 


and 


SDP 


where  x denotes  the  value  of  x evaluated  at  the  saddle  point  a =9. 
S s 
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(92) 


-k2  fdae"jkDcos(a-9)sin2avl(x)]  . (93) 


(94) 


dae-jkDcos(a-9)sin2aV1(x)~V1(xs)sin26  e‘jkD  (95) 


— - 


In  view  of  the  relation  (35)  and  (37)  one  finds 


f 2 w cos2/30  Me'ikD= 


e J 
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„ , v . 2q  / 2TTi  -ikD  . „„,3  , x sin20  e'^0 

Substituting  (96)  and  (97)  into  (93),  one  has  approximately 


jM 

TT»C,  /\  ! ..CD 

Hcp  (£*£  > 2 

^ ' / 4it  mu  a 
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2nka7TLvo(xs)£Eir-j 
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sin29  e -jkD  "j 


+ J2nk-avl(*a)^=— yj  (S 

cos  O (kD)  I 

Rigorously,  the  asymptotic  evaluation  of  the  integral  in  (94)  should 
be  carried  out  to  next  order. 

Since  for  x < 1,  v (x  ) is  slowly  varying  (see  Figure  II-4),  one  may 
u s 

employ  the  approximation 

d2  T , x e'jkD1  . . 92  e'jkD 

k2v0(%)(3in2e-3in29ife  + cos29lt”i)£ki5n'  <■ 

where  after  differentiation  terms  o( — ^r)  have  geen  neglected.  Meanwhile 

/ i > (kD  )J/ 

the  next  term  in  (95)is  0( 575-]  in  the  deep  shadow,  so  that  in  view  of  (37) 

(kD)3/2/  j 

the  corresponding  near  zone  expression  will  be  o(^^j3  J and  is  neglected. 

Substituting  (99)  into  (98),  one  finds  the  approximate  field  expression 


Hjp  (r.r')|  « Mj,  ^r[“j  (sin20  - sin29  |jj  ) 


p=p  =a 


sin^e  1 x^  e'jkD 

.2a  IcD  V1  Xs  J kD 


The  soft  contribution  is  obtained  directly  from  (78).  Combining 


(78)  with  (100)  one  finds  the  "Full  Formula" 


H^(£.r ')  | »-j  Mr^(sin20  - sin28-L  + cos20  e 
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cos  o (kD) 

The  properties  of  this  Full  Formula  are  as  follows: 

a)  it  matches  to  of — W/the  asymptotic  formula  (79)  for  0=0,  i.e. 

(kDr 

in  the  circumferential  (H-plane), 

b)  it  is  finite  on  0 = and 

c / i \ 

c)  as  D ■*  0 H,  tends  to  the  planar  result  to  O ( =■  j . . 

' /I 


(kD) 


This  is  seen  as  follows: 


w(xs)  = v0(Xs}  + 


sin29 Vj(xg)  -u0(xg) 

re 

COS  0 


f2(l+  3Si„26)  - + sin29)  cos29+  O (gf 


(101) 


(102) 


Consequently  w(xg)  “*  0 as  D - 0 and  (101)  reduces  to  the  planar  result  for 
H^,  where  C denotes  the  Cartesian  coordinate  in  the  direction  of  the  magnetic 
dipole  on  a ground  plane  and  9 is  the  angle  of  the  radius  vector  with  the  £ 
direction. 

Following  the  same  line  of  reasoning  one  has  from  (12)  and  in  the 
first  quadrant 


H^r.r') 


j M, 
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After  certain  amount  of  algebra  one  finds 


H=(r,r') 


0=  P =a 


M Y d2  -jl 

v*.1 57v  -n 


where  the  integrals  were  evaluated  to  lowest  order  in  t-=t  and  vn(x  ) was 


assumed  slowly  varying.  Hence  approximately 


H'(r.r') 


• w k2Y  .0  af,  3j  ''  , v e'jkD 

J Irr  9in  9 cos  9 ; 1 - vQ(xg)  _w_ 


0=P  =a 


b.  Variant  of  the  ' Full  Formula' 


This  formula  is  obtained  from  (101)  by  retaining,  in  addition  to  the 


soft  polarization  contribution  Uq(x),  only  hard  terms  to  lowest  order  in 
in  line  with  the  rigorously  asymptotic  procedure  used  in  (74).  Examining  in 


the  deep  shadow  the  hard  polarization  terms  Vq(x)  and  v^(x)  one  observes,  in 


view  of  (35)  and  (37),  that  for  a given  distance  D 


vi(It)  0r  i 

CWJ7> 


since  from  (62)  and  (63)  one  has  x = k^  = — ~J 

2£g 


Thus,  for 


consistency,  the  term  involving  v^(xg)  should  be  discarded.  As  a result,  in 
the  paraxial  region,  for  xg  <0.8  or  0.  6,  it  is  necessary  to  neglect  the  soft 
polarization  term  Uq(x),  which  becomes  unbounded  for  9 ^ • 


Thus,  the  variant  formula  becomes 
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It  should  be  noted  that  in  the  0=0°  (H-plane)  both  the  ’’Full  Formula" 
and  its  variant  coincide  with  the  rigorously  asymptotic  result  of  (79). 
c.  Axial  Magnetic  Current  Element 
From  (8),  (11),  (16)  and  (24)  one  has 
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In  a similar  manner,  retaining  in  the  integral  of  (10 9)  only  the  leading 
term,  one  has 


v f ? f"  "jkDi 
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(1 


0=0  =a 


Again,  assuming  that  vn(x  ) is  slowly  varying  for  x < 1,  one  finds 

US  s 
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In  view  of  the  identical  form  of  (25)  and  (12)  and  because 


3z  dcp 
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dz  ' dro ' 


C" , one  has,  via  (105) 
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• jkD 
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The  range  of  9 in  (101)  and  (105)  and  (11-112)  is  0 < 9 < 2tt,  i.e. 
over  all  four  quadrants.  This  feature  automatically  accounts  for  the  sign 
changes  from  quadrant  to  quadrant,  as  required  by  (80)  and  (91)  . 

d.  Modified  Planar  Vector  Formula  for  Arbitrary  Orientation 
of  Current  Source 

The  set  of  approximate  relations  (101),  (105),  (111)  and  (112)  can  be 
cast  in  a vector  form  similar  to  that  for  the  planar  case  (see  Appendix  D). 
In  Figures  II- 10(a)  and  (b),  the  unit  vectors  t and  t*  are  tangent  vectors, 

6 and  6 binormals  to  a typical  surface  ray  and  s denotes  the  source 


(a)  Circumferential  M_  (b)  Axial  M 

Figure  11-10.  Surface  Rays  on  a Developed  Cylinder 


direction.  The  prime  indicates  source  coordinates  while  the  unprimed 

quantities  refer  to  the  observation  point.  The  angle  9 is  given  to  tan’  — 

a(ffl-a) ') 

One  now  obtains  for  a circumferential  magnetic  current  element  (Figure  II- 10a); 

A A 

i.  e.  s = in 
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Therefore, 
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The  second  term  in  the  bracket  in  (105)  may  also  be  written 
sin  9 cos  9 + 2 cos  9 sin  9,  Therefore,  one  has,  via  (115), 
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Making  use  of  the  unit  dyadic  cntn  + z z = 1,  one  finds 
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For  an  axial  M,  (see  Figure  II- 10b) 
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cos9sin0=-s*b  6 • sp  = t • cp  s • t . 

Hence,  in  a similar  manner  as  in  (118),  one  has 
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Combining  both  slot  orientations,  one  finally  has 
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where  M denotes  an  arbitrarily  oriented  magnetic  point  source.  Note, 
this  formula  is  identical  with  the  planar  formula  (D8),  save  for  the  curvature 
terms  vQ(xg)  and  w(xs)* 

The  analysis  of  the  circular  cylinder  problem  was  carried  out  in 
terms  of  the  Cp  and  z coordinates  which  simultaneously  constitute  the  direc- 
tions of  principal  curvature  of  the  surface.  In  its  generalization  to  a conical 
surface  V - 9 and  z - r,  » and  r again  being  the  principal  curvature  directions 

on  a cone  (see  Chapter  IV). 

For  a more  general  surface  of  revolution,  the  tangential  M will  have 
to  be  decomposed  prior  to  further  treatment  into  components  along  the  two 
principal  directions. 

6.  Summary  of  Formulae 

a.  Strictly  Asymptotic  Expressions  to  Lowest  Order  in  ^ and 
to  Second  Order  in^: 

A.  Circumferential  Magnetic  Current  Element 


42 


H.^(r,r')|  * - j sin20vn(x_)  e 


r 2^ 


s 


kD 


0=0  =a 


. k2Y.'  2 a 61  . 2fl  1 sin26N  , v e'jkD 

- 2TT  " 24  9 0 3 (kD)2 


+ M 


k Y 1 


e-jkD 


I7u0(xs)— ITZ 

cos  o (kD) 


,2  v -jkD 

Hz  (i*  X #)  | ~ J M-TfTsin0  cos8v0(xs)^kD- 

0=0 '=a 


w k2Y7  23  . o a . 5 sin09  e 

- M^j  ~2rT \ ~ IT  sin  ^ Cos9-'-9'^rB/v0(xs)- 


-jkD 

T 


(kD) 


B.  Axial  Magnetic  Current  Element 


i 2.,  , -jkD 

r*a  . /.  | . . , k Y 20  ,.eJ 

Hz(r,r  )|  --JM^cos  9 W Td” 


0=0  =a 


, 2V  -jkD 

H*(lV)|  - j Mz^sin9  cos  9 vQ  (xg) 

0=0 '=a 


- M 


kY{  23 


Validity  of  these  formulae  is  restricted  to  9 not  too  close  to  11 /2  , 
b.  The  Modified  Planar  Formula 

This  formula  is  not  strictly  asymptotic.  It  applies  to  an  arbitrary 
M,  and  is  written  in  a vector  form. 


(124) 


(125) 


(126) 


(127) 
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n^L)\  ~ * j 7tt m*  l6'6(1  - is)  + v*lfc>o(x8)i 

-jkD 


-jkD 

EE- 


P=P  =a 


- M-  $'  ^rT  w(xg)  - 


where 


W(X3)  = W + 


(kD) 


sin  0 v (x  ) - u (x  ) 
Is  os 

2fl 
COS  o 


(128) 


(129) 


c.  The  "Full  Formula" 

The  expression  for  excited  by  a circumferential  magnetic  current 
element  derived  in  (101)  or  obtainable  from  (128)  - (129)  is  termed  here  the 
"Full  Formula. " 


2 

K ti')  I - - i M;p  irr(sin2e  - sin2e  + cos20  |^;v0(xg) 

.-jkD 


■jkD 


kD 


p=c 


w k2Y  , , e 

* Nl...  1 W X ) 

CP  ZTT  3 


(kD) 


(130) 


where  w(xg)  is  shown  in  (129). 


d.  The  Variant  of  the  "Full  Formula1 

2. 


«£  (r,  r #)  | ~ - j ^/(3in28  - sin26  ^ + cos20  Vxg> 


.-jkD 


0=0 


. w k2 Y 1 , , e"jkD 

+ Mcplrr^T7u0(xs)7^T 


cos  0 


(kD) 


for  x > 0.  6 
s 


• ™ k“Y / . 2Q  . 2fl  2i  . 2fl  j \ . . e'jkD 

--■»Mq>-Srv*ln  9'Sin  9kD+cos  9kD^v0(xs)-TE~' 


for  x <0.6 
s 


(132) 


1/3 


where  x 


kD  \ a 

s ~F’  £g  * ’vTZ*  j and  pg  = 27  <see  62  and  63>- 

2f  6 s cos  9 

g 

The  curvature  terms  vo(Xg),  v^x^and  uQ(xg)are  listed  in  (35-43). 
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III.  Numerical  Results 


This  chapter  is  devoted  to  numerical  tests  of  validity  of  surface 
ray  methods  for  mutual  coupling  in  conformal  arrays.  Specifically,  these 
methods  are  applied  to  mutual  coupling  between  slots  on  conducting  circular 
cylindrical  surfaces.  Once  the  limits  of  validity  of  the  surface  ray  tech- 
niques have  been  established,  the  asymptotic  results  of  the  cylindrical  prob- 


lem can  easily  be  extended  in  a GTD  manner  to  conical  surfaces.  This  is 
described  in  Chapter  IV. 

A literature  search  has  disclosed  the  availability  of  numerical  and 
experimental  results  for  mutual  admittance  between  "single  mode"  rectangu- 
lar  slots  in  a conducting  circular  cylinder  [5],  [9].  These  were  obtained  via 
harmonic  series  in  X and  a numerically  evaluated  Fourier  integral  with  respect 
to  the  axial  wave-number  spectrum.  To  enforce  convergence,  a small  loss 
was  assumed  in  the  surrounding  medium.  The  authors  were  helpful  in  sup- 
plying us  with  some  of  their  numerical  data  in  a computer  print-out  form. 

They  also  made  available  their  computer  program  to  Drs.  Kummer  and 
Villeneuve  of  Hughes  Aircraft  Co.,  Culver  City,  California,  who  in  turn 
furnished  us  with  additional  numerical  harmonic  series  results. 

The  surface  ray  computations  of  mutual  admittances  between  "single 
mode"  rectangular  slots  in  an  infinite  conducting  circular  cylinder  were 


mm 


normalized  so  that 


rr 

j .■ 


—2 


dS  = 


Aperture  2 


rr 

v O 

Aperture  1 


dS  = 1 


(2) 


and  H21  is  the  surface  magnetic  field  intensity  due  to  aperture  1 on  the 
unperforated  surface  in  the  location  of  aperture  2.  The  surface  magnetic 
current  density  is 


M = E x n dS  (3) 

where  n is  the  local  outward  unit  normal  of  the  cylindrical  surface  and  dS 

denotes  an  infinitesimal  area  in  the  slot.  In  the  case  of  circumferential  slots, 

it  follows  from  (1)  and  (3)  that  only  needs  to  be  considered,  while  for 

axial  slots  only  H is  required. 
z 

To  compare  to  Stewart-Golden1  s data,  the  slot  dimensions  were 
chosen  those  of  a standard  X-band  0.9"  x 0.4"  waveguide. 

a.  Circumferential  Slots 

Two  formulae  for  were  used  in  obtaining  the  numerical  values 
for  Y^*  the  "Full  Formula"  (11-101)  and  its  variant,  (11-107,  108)  as 
described  in  Section  5 of  Chapter  II. 

In  the  following,  Zq  and  i denote  the  axial  and  the  angular 
center-to-center  separation  of  two  circumferential  slots  (see  Fig.  Ill  - 1 . 


For  further  comparison,  curves  based  on  Hwang  and  Kouyoumjians 
GTD  formulation  [10]are  also  included.  These  formulae  consist  solely  of 
the  dominant  hard  and  soft  terms  of  (11-123).  For  a circumferential 
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Figure  III  — 1 . (a)  Two  Circumferential  Slots  on  a Cylinder  of  Radius  a = R, 

(b)  Slot  Geometry, 

(c)  Partitioning  of  Slots  for  Numerical  Integration, 
magnetic  current  element,  in  our  notation,  these  formulae  are 

,2  _ - jkD 

P=  P '=a 

k^Y  1 -jkD 

+ M — Uj-  un(x  ) * 


and  for  an  axial  magnetic  current  element, 

, 2V  ? _ - jkD 

~ -jMz'2^rcos  9v0(xs)"TD_  • (5) 

p=  p'=a 

We  shall  consider  the  mutual  admittance  between  circumferential 

0.  9"  x 0.4"  single  mode  slots  in  conducting  cylinder  with  a = R = 1.  991"  at 
f = 9 GHz  (ka  » 9.  5)  and  subsequently  with  a = 3.  777"  at  f = 0.  975  GHz 

(ka  « 19.  6).  In  each  case  the  data  are  normalized  to  the  magnetic  of  the  self 
admittance  |yu1,  where  Yj  j / Yg  = 0.  8588  + j 0.  3834  for  a = 1 . 991"  at 
9 GHz  and  Y,  , / Y = 0.  8440  + j 0.  4040  at  9.  75  GHz, for  a = 3.  777",  Y being 

1 i g g 

the  TE1q  modal  admittance  of  the  rectangular  feed  waveguide  with  0.9"  x 0.4" 

1.  D.  The  computed  points  are  in  10°  steps  of  angular  separation  $ for 
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a = 1.  991"  and  in  10.  1333°  for  a = 3.  777". 


1 


The  procedure  of  numerical  integration  of  (1)  for  circumferential 
slots  was  based  on  a subdivision  of  each  slot  into  14x2  = 28  elemental 
areas  AS  of  equal  size  (Figure  III-lc).  An  increase  of  the  number  of  sub- 
divisions to  14x4  = 56  reduced  the  difference  between  S-G  and  PINY  by 
0.  3 DB  for  the  near  slots  ($q  < 50°)  but  did  not  appreciably  affect  the 
numerical  results  for  the  far  apart  slots. 

We  shall  begin  with  the  H-plane  coupling  (Zq  = 0)  and  increase  Zq 
in  steps.  Figure  (III- 2 ) shows  jY^/Y^jl  calculated  by  the  Full  Formula 
(PINY)  for  ZQ  = 0,  in  comparison  with  the  harmonic  series  data  (S-G), 
and  GTD  . The  PINY  curves  show  good  agreement  with  S-G, 

but  the  GTD  curve  markedly  deviates  from  S-G  and  shows  a different  slope. 
As  discussed  in  Chapter  I,  the  reason  for  this  discrepancy  is  that  the  lowest 

order  in  hard  contribution  vanishes  for  9=0  and  the  behavior  of  H„  is 
kD  cp 

dominated  not  by  the  soft  contribution,  but  rather  by  the  second  order  in 

hard  contribution,  the  exponential  attenuation  rate  of  which, is  identical 

with  that  of  the  lowest  order  hard  contribution,  the  attenuation  constant  of 

the  dominant  hard  surface  ray  being  that  of  the  dominant  soft  surface 

ray  for  the  same  geodesic.  Figure  III- 3 shows  the  effect  of  neglecting  the 

Vj(x)  term  in  (II- 101)  and  in  addition  discarding  the  soft  contribution  Up(x) 

for  x < 0.  3 (actually  x = 0.  6 would  be  more  appropriate).  In  the  H-plane 

sin20 

the  Vj(x)  term  has  a negligible  effect,  because  its  pattern  has  a 

cos  9 

null  for  9=o. 

Table  I gives  the  phase  comparison.  It  is  seen  that  the  "Full 
Formula"  (F.  F. ) yields,  except  for  = 11C°,  a maximum  phase  difference 
of  » 13°  from  S-G. 
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1 


i 

/ 


! 


Figure  III-4  shows 


12 


11 


for  the  case  of  Zq  = 0.5".  The  departure 
of  the  F.  F.  from  the  S-G  result  increases  and  the  departure  of  GTD  decreases. 
Table  II  summarizes  the  phase  comparison  and  it  is  seen  that  the  correlation 
of  the  F.  F.  to  S-G  is  very  good,  next  comes  the  variant  formula  and  finally 


GTD. 

Figure  III- 5 shows  the  Zq  = 1"  case.  The  F.  F.  curve  is  now  above 
S-G  and  so  is  the  variant  formula  which  shows  now  a better  amplitude  cor- 
relation than  F.  F.  with  the  harmonic  series  result.  The  phase  comparison 
is  shown  in  Table  III.  Again  best  phase  correlation  is  obtained  via  F.  F. 

Figure  III-6  and  III- 7 and  Table  IV  present  the  case  of  Zq  = 2".  The 
F.  F.  is  about  2DB  higher  than  S-G,  while  GTD  correlates  better  for  $q>30°, 
but  tends  to  diverge  for  small  values.  On  the  other  hand  the  variant 
formula  result,  to  the  accuracy  of  the  graph,  is  indistinguishable  from  S-G 
except  towards  the  E-plane  where  it  falls  below  S-G  by  about  1 DB.  Similar 
comments  apply  to  the  case  Zq  = 3",  represented  Figure  III- 8 and  III- 9 and 
Table  V.  It  is  seen  that  in  Figure  III-8  the  difference  between  the  F.  F.  and 
the  harmonic  series  result  is  -ZDB.  On  the  other  hand  the  variant  formula 
and  GTD  track  S-G  quite  well. 

Figures  III- 10  and  III- 11  and  Table  VI  summarize  the  comparison  for 
the  case  of  Z = 4".  There  is  about  1.7DB  difference  in  the  E-plane  for 
either  F.F.  or  the  variant  formula.  The  F.F.  phases  track  very  well,  next 
best  is  the  variant  formula. 

Figure  III- 12  and  Table  VII  present  the  results  for  Z = 8".  The 
maximum  deviation  of  the  F.F.  result  from  that  of  harmonic  series  is  about 
3 DB.  However,  it  is  to  be  noted  that  in  view  of  the  assumed  loss,  i.e. 
e = Sq(1  - j 0.  0025)  the  actual  discrepancy  should  be  reduced  to  no  more  than 


■ 


2 DB. 
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As  a general  comment,  the  GTD  formulae  show  good  agreement 

in  amplitude  with  the  harmonic  series  for  Z > 2"  and  $ > 30°.  The  variant 

o o 

formula  matches  very  well  to  S-G  in  amplitude  over  the  entire  range  of  Z 

o 

including  the  H-plane,  but  falls  below  the  harmonic  series  result  in  the 
E -plane.  Also  the  phase  of  the  variant  formula  tracks  better  that  of  S-G 
than  the  GTD  result.  The  phase  of  the  F.F.  tracks  best  of  all.  The  dis- 

Y12  I 

crepancy  in  yz — is  increasing  with  Z , but  seems  to  level  off  for  large 
1 1 1 1 ° 

values  of  Zq,  in  view  of  the  loss  assumed  in  S-G  program. 

Figures  III- 13  to  III- 17  and  Tables  VIII  to  XII  show  the  amplitude 
y12 

and  phase  of  -n — for  a = 3.  Ill" , f = 9.  75  GHz  and  0.  9"  x 0.  4"  apertures 
*11 

for  different  axial  slot  separations  ZQ.  In  this  case,  the  harmonic  series 
computer  program  employed  a loss  tangent  of  0.  003.  As  expected,  for  this 
increased  value  of  ka<»19.  6,  the  correlation  between  the  harmonic  series 
and  the  asymptotic  surface  ray  results  is  improved  in  amplitude, 
b.  Axial  Slots 

Since  the  expressions  (11-90)  are  finite  as  6 they  may  be 

adequate  for  the  description  of  the  field  in  the  neighborhood  of  a slot  uni- 
formly in  all  directions  and  thus,  suitable  for  the  determination  of  the 
mutual  admittance.  One  observes  that  for  axial  slots  only  the  hard  polar- 
ization is  excited,  so  that  the  soft  transition  function  uq(x)  does  not  appear 
in  (11-90). 


Samples  of  numerical  results,  obtained  via  harmonic  series,  for 
mutual  coupling  between  single  mode  axial  slots  on  a circular  cylinder  can 
be  found  in  [9].  These  data  are  in  the  form  of  isolation  20  logjQ  Is^  1 
between  two  open-ended  rectangular  waveguides  feeding  the  slot  apertures, 
where  S^  is  the  scattering  (coupling  coefficient  between  the  waveguides). 
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4>  DEGREES 
o 

Figure  III-3.  Mutual  Admittance  of  Circumferential  Slots 

on  a Conduetina  Cylinder  a =__l._991"i  £0=.Q . . 


TABLE  I Phase  Comparison,  Arg.  Y 

Z = 0,  R=  1 . 991  ’ , f=9  GHz 
o 


No.  Vj 
DEGREES 

Full  Formula 
DEGREES 

GTD 

DEGREES 

S-G 

DEGREES 

o 

3 

- 

- 

- 

- 

0 

- 

- 

- 

- 
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- 

- 

- 

- 
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-8.  21 
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-103. 25 

-153.  17 
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153. 70 


152.  82 


92.71 
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TABLE  II 


Phase  Comparison  Arg.  Y 
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Z =.  5",  R=  1.  991",  f=9  GHz 
o 


No  V1 
DEGREES 

Full  Formula 
DEGREES 

GTD 

DEGREES 

S-G 

DEGREES 

gP 

0 

-81.  84 

-70.  85 

-72.  38 

10 

-85.  57 

-77. 63 
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20 

-99. 46 

-107. 29 
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30 

-151.  36 
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-144. 44 
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40 
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50 
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8.  77 

2.  51 

7.  45 

60 

-81.  43 
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70 

176.  65 

168. 89 
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168. 20 

80 
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-28. 18 
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-41.  51 
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-178. 93 

-142.  26 
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127.49 
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73.  03 

131/22 

TABLE  III 

Phase  Comparison, 

Arg.  Yu 

Z =1”  , R = l.  991 
o 

",  f=9  GHz 

No  Vj 

Full  Formula 

GTD 

S-G 

DEGREES 

DEGREES  DEGREES 

DEGREES 

o 

9 

0 

152.  62 

152. 99 

- 

154. 67 

10 

151.  24 

142. 84 

- 

144. 69 

20 

134. 53 

111. 19 

- 

113. 51 

30 

86.  51 

57.  95 

101. 53 

61.  20 

40 

9.  01 

-12.  30 

13.  62 

o 

O' 

r- 

i 

50 

-75.  81 

-93. 86 

-77. 03 

-88.  37 

60 

-165.  88 

177. 49 

-172.  19 

-176. 18 

70 

100. 36 

84.  48 

89.  47 

91.  39 

80 

4.  22 

-11.  27 

-10. 91 

-3.  58 

90 

-93.40 

-108. 69 

-112. 58 

-100.  29 

100 

168.  07 

152. 88 

144.  99 

161.  02 

110 

69.  00 

53.  87 

42.  12 

61.  96 

MUTUAL  ADMITTANCE  OF  CIRCUMFERENTIAL  SLOTS 
ON  A CONDUCTING  CYLINDER 


ELIMINATE  v,  (x)  AND  LEAVE  OUT 
u0(x  ) FOR  x < 0.8 
STEWART-  GOLDEN 


Figure  III  - T . Mutual  Admittance  of  Circumferential  Slot 
on  a Conducting  Cylinder,  a = J.9S-L ?,Z-s  - 


TABLE  IV 


Phase  Comparison,  Arg.  Y 
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Z =2' 
o 


R=l.  991",  f=9  GHz 


No  Vl 

Full  Formula 

GTD 

S-G 

DEGREES 

DEGREES 

DEGREES 

DEGREES 

o 

CP 

0 

-113. 06 

-118.21 

- 

-117.  01 

10 

-120.75 

-126. 10 

(-10.  33) 

-124. 73 

20 

-143.46 

-149.45 
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30 

-178. 14 
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175. 34 

40 
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9.  25 
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-86. 36 
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MUTUAL  ADMITTANCE  OF  CIRCUMFERENTIAL  SLOTS 
ON  A CONDUCTING  CYLINDER 


STEWART -GOLDEN 


Mutual  Admittance  of  Circumferential  Slots 
on  a Conductine  Cylinder-  a = I.  991 M,  Z = ’ 


TABLE  V Phase  Comparison,  Arg.  Y 


7 - V ' 
o 3 ' 

No  V 

DEGREES 

R=l.  991",  f=9  GHz 
Full  Formula 
DEGREES 

GTD 

DEGREES 

S-G 

DEGREES 

o 

% 

. 

0 

-22.  89 

-32.  33 

- 

-30.  90 

10 

-28. 67 

-38.  19 

(75.  84) 

-36.  72 

20 

-45.  76 

-55.  67 

(72.  46) 

-54.  11 

30 

-71.  91 

-84.  27 

-22.  20 

-82.  33 

40 
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-123.  14 

-83.  23 

-120. 52 

50 

-147. 55 

-171.  15 

-140. 36 

-167.  78 

60 

153. 96 

132.  82 

+158. 90 

137. 00 

70 

88.94 

69.  90 

93.  13 

74.  88 

80 

19.01 

1.  15 

22.  42 

6.  78 

90 

-55. 24 

-72.  46 

-52.  59 

-66. 24 

100 
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-150. 12 

-131. 53 

-143.  17 

110 

145. 59 

129. 69 

146. 75 

137.  29 

MUTUAL  ADMITTANCE  OF  CIRCUMFERENTIAL  SLOTS 
ON  A CONDUCTING  CYLINDER 


GTD  +16.6 +3dB 


ELIMINATE  v,  (x)  AND  LEAVE  OUT 
u0(x ) FOR  x <0.8 
STEWART- GOLDEN 


Mutual  Admittance  of  Circumferential  Slots 
on  a Conducting  Cylinder,  a=l.Qal".Z  =<l 


Figure  III- 11 


TABLE  VI 

No  V1 
DEGREES 

Phase  Comparison, 

Z =4",  R= 1. 991 
o 

Full  Formula 
DEGREES 

Arg.  Y12 
",  f=9 GHz 

GTD 

DEGREES 

S-G 

DEGREE 

o 

CP 

• 

0 

65.  14 

52.98 

- 

54.  11 

10 

60.  54 

48.  33 

(161.  79) 

49.  54 

20 

46.  83 

34.  45 

(158.99) 

35.  73 

30 

24.  31 

11.  57 

81.  11 

13.  08 

40 

-6.  57 

-19.81 

25.  54 

-17.  72 

50 

-45. 21 

-58.  98 

-24.  85 

-56.  21 

60 

-90. 97 

-105. 25 

-76. 76 

-101. 83 

70 

-143.  16 

-157. 94 

-132. 75 

-154. 02 

80 

159.01 

143. 71 

166. 84 

147.  94 

90 

95.  98 

80.  36 

101.  94 

84.  99 

100 

28.  81 

12.  66 

33.  31 

18.  78 

110 

-42.  76 

-58.  53 

-39.  34 

-50.  84 

TA  BLE  VII 

Phase  Comparison,  Arg.  Y. 

Z =8",  R=l. 
o 

991",  f=9GHz 

No  Vx 

Full  Formula 

GTD 

S-G 

DEGREES 

DEGREES 

DEGREES 

DEGREES 

D° 

0 

51.61 

33.42 

- 

36.  67 

10 

56.  82 

30.  84 

- 

30.  09 

20 

71.02 

23.  10 

- 

22.  65 

30 

78.  17 

10.  29 

98.  60 

10.  18 

40 

53.  09 

-7.66 

56.43 

-7.  33 

50 

15.  30 

-30. 32 

17.  55 

-29.  57 

60 

-21.  60 

-57. 52 

-19.  13 

-56.  33 

70 

-58.  87 

-89.05 

-56. 35 

-87. 29 

80 

-98.  08 

-124. 69 

-95. 58 

-122.47 

90 

-139.91 

-164. 18 

-137.45 

-161.  78 

100 

175. 38 

152.71 

177.77 

155. 71 

110 

127. 73 

10.  19 

130. 05 

70.  71 

1 


TABLE  VIII  Phase  Comparison,  Arg. 

Z =0,  R=  3.  777”,  f=9,  75GHz 
o 


No  V1 

Full  Formula 

GTD 

S-G 

DEGREES 

DEGREES 

DEGREES 

DEGREES 

CfP 

V 

- 

_ 

0 

- 

- 

- 

- 

10.  13 

- 

- 

- 

- 

20.  27 

178. 24 

177. 92 

176. 11 

166.  5 

30.40 

-45.  28 

-45. 63 

53.  53 

-54. 65 

40.  53 

97.  81 

97.  55 

82.  59 

89.  20 

50.  67 

-116.  29 

-116. 52 

-140.40 

-124.  89 

60.  80 

31.46 

31.  26 

-3.  36 

27.  71 

70.  93 

180. 73 

180. 55 

133. 64 

167. 26 

81.  07 

-28.  80 

-28.  96 

-89. 23 

-29.  50 

91.  20 

-122.44 

122. 29 

48.  11 

105. 03 

101.  33 

-85.  98 

-86.  11 

-173. 96 

-77.49 

111.47 

65.  68 

65.  46 

-35. 23 

9.97 

! 


J 
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TABLE  IX  Phase  Comparison,  Arg. 

Z =.455",  R=  3.  Ill",  f=9.  75GHz 
o 


No  Vx 
DEGREES 

Full  Formula 
DEGREES 

GTD 

DEGREES 

S-G 

DEGREES 

5° 

0 

-81.  78 

-64.43 

-66.  36 

10.  13 

-111.  94 

-98. 63 

- 

-98.  84 

20.  27 

118. 33 

109. 17 

112. 06 

105. 25 

30.  40 

-83.  20 

-87.79 

-100. 20 

-92.  43 

40.  53 

69.35 

65.99 

44.  90 

61.  76 

50.  67 

-139. 68 

-142. 50 

-171.  76 

-146. 67 

60.  80 

11.  22 

8.  70 

-29. 22 

3.  66 

70.  93 

162.64 

160. 34 

113. 27 

158. 34 

81.  07 

-45. 27 

-47.  37 

-103.  72 

-57. 03 

91.  20 

107. 33 

105.  42 

40.  24 

109. 48 

101.  33 

-99. 88 

-101. 58 

-174. 55 

-116. 88 

111.47 

52.  80 

51.  28 

-27.  91 

65.  49 
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TABLE  X 


Phase  Comparison, 

Arg  Y12 

7 - 2n 

o'* 

R=  3.  777",  f 

= 9.75  GHz 

o 

No  V1 
DEGREES 

Full  Formula 
DEGREES 

GTD 

DEGREES 

cp 

0 

-157.72 

-159.06 

- 

10.  13 

178.45 

171.  24 

- 

20.  27 

109.  69 

86.  80 

117.89 

30.40 

-25. 78 

-40. 14 

-21.49 

40.  53 

176. 63 

164.75 

177.88 

50.67 

3.  11 

-7.  87 

1.  87 

60.  80 

179. 35 

168.71 

175. 27 

70.  93 

-10.92 

-21.46 

-17.40 

81.07 

154.  59 

144. 02 

145. 82 

91.20 

-42.70 

-53. 37 

-53.70 

101.  33 

118.  09 

107. 29 

104. 92 

111.47 

-82. 52 

-93.46 

-97.81 

S-G 

DEGREES 
-158. 07 
173.22 
91.24 
-32.46 
176. 39 
1.  58 
-171.  96 
0.  29 
171.41 
-27. 29 
145. 33 
-69.  45 
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r 


Phase  Comparison, 

Arg.  Y12 

7 =4” 
o * 

R=3.  777",  f=9. 

75  GHz 

No  Vl 
DEGREES 

Full  Formula 
DEGREES 

GTD 

DEGREES 

S-G 

DEGREES 

0 

-25.  69 

-32. 26 

- 

-31. 74 

10.  13 

-35.  18 

-49. 52 

- 

-48. 15 

20.  27 

-67. 04 

-100.03 

-58. 55 

-97. 04 

30.  40 

-160.  52 

-179.  47 

-156. 31 

-175.  26 

40.  53 

88.  30 

73.  46 

91.  82 

81.  67 

50.  67 

-40. 36 

-53. 06 

-37. 61 

-41.  78 

60.  80 

175.  89 

164. 21 

177.  86 

179.  08 

70.  93 

20.  10 

8.  94 

21.  29 

27.  21 

81.07 

-145.  11 

-155.  98 

-144. 67 

-134. 39 

91.  20 

42.  35 

31.  61 

42.  05 

57.  36 

101.  33 

-135.  96 

-146.  62 

-136.97 

-118.  90 

111.47 

41.  15 

30.  51 

39.45 

64.  32 

' 

i 
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TABLE  XII  Phase  Comparison,  Arg. 

Z =8",  R = 3.  777",  f=9.  75GHz 
o 


No  VL 
DEGREES 

Full  Formula 
DEGREES 

GTD 

DEGREES 

S-G 

DEGREES 

o 

& o 

-131.01 

-142. 09 

. 

-140. 77 

10.  13 

-133. 94 

-151.40 

- 

-150. 35 

20.  27 

-134. 94 

-179.04 

-119. 36 

-177. 69 

30.  40 

167. 65 

135. 71 

169.93 

138. 36 

40.  53 

95.  66 

73.  59 

98.  00 

78.  16 

50.  67 

13.91 

-3.74 

16.  15 

2.  89 

60.  80 

-79. 61 

-94. 97 

-77. 51 

-85.  86 

70.93 

175.  36 

161.  32 

177. 30 

173. 67 

81.07 

59.  72 

46.  51 

61.48 

61.  87 

91.  20 

-65.48 

-78. 13 

-63.  91 

-60. 59 

101. 33 

160. 79 

148. 54 

162. 17 

170. 44 

111.  47 

19.47 

7.  53 

20.  66 

32.  90 

1 
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ter 


Figure  III- 18.  Axial  Slots  on  a Cylinder 
To  arrive  at  a comparison,  it  was  necessary  to  convert  the  values 

obtained  via  (1)  and  (3)  to  via 

-2  Y Y12 

S12  = (Yg+  Vn  + V6u)(Yg+Vn.  V12)  " 

where,  again  Y is  the  characteristic  admittance  of  the  dominant  feed- 
g 

waveguide  mode  and  Y^  the  slot  self-admittance.  To  save  computer  time 
the  number  of  subdivisions  for  the  axial  slot  was  chosen  10x2  = 20,  and 
the  tests  with  14  x 4 = 56  subdivisions  changed  the  result  by  0.  1 DB  even 
for  the  nearby  slots. 

Figure  III- 1 9 and  III-20  show  the  comparison  between  the  harmonic 
series  and  the  surface  ray  results  of  E-plane  coupling  (ZQ  = 0)  and  Zq  = 1.5" 
at  9 GHz,  for  ka  « 9.  5.  A very  good  agreement  is  found  everywhere  within 
the  angular  separation  of  20°  to  150°  for  the  case  of  Zq  = 0,  and  0°  to 
150°  for  Z = 1.5".  The  GTD  formula  is  completely  inadequate  in  the 
H-plane  (axial  direction). 

Y11 

From  [9],  = 0.  748  + j 0.  434. 
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IV.  GTD  Generalization  to  Conical  Surfaces 

a.  GTD  Generalization  of  the  Curvature  Factors  to  Surfaces 
with  Variable  Curvature 

The  generalization  will  be  discussed  for  vQ (x),  the  other  curvature 
terms  are  treated  in  a similar  fashion. 

For  a circular  cylinder  one  has  from  (11-36) 


. TT  ® -jXt 

-J  4 V * P 


/ » J 4 r—  r—  \ e P 
Vq(x)  = e v rr  ,/x  ^ — — , in  deep  shadow, 

p=l  *p 


1 x3/2  + {^x3  + - x9/2+. . . , in  the  near  zone,  (lb) 

VT  512 


where 


kD  , (^z) 

x = rpr  * fg  = hr  J * pg  = — 27 

2 fg  6 cos  o 

and  Pg  is  the  principal  radius  of  curvature  of  the  ray  trajectory.  On  a 

circular  cylinder  is  constant  along  a ray.  On  a cone  however,  or  other 

surfaces  of  variable  curvature  P = P (s),  varies  along  the  ray,  where  s 

S 8 

is  distance  parameter  along  a ray. 

The  generalization  of  v^(x)  to  a conical  surface  proceeds  as  follows. 
Firstly  x in  the  exponent  is  replaced  by  x such  that  for  an  infinitesimal 
distsnce  along  a ray 


d x(s ) = 


and  consequently 


— |*  k ds 

x = I — x — 


2f  (s) 
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i 


In  (4)  O'  denotes  the  source  point,  Q the  observation  point  and  the  ray- 
trajectory  passes  along  a geodesic  path  from  Q ' to  Q. 

Secondly,  in  (la),,/x  which  represents  the  product  of  the  launching 
and  attachment  coefficients  [l],is  symmetrized  in  order  to  preserve  recipro- 
city with  respect  to  an  interchange  of  the  source  point  (Q  7)  and  observation 
point  (Q).  Thus, 


~J — , 

V 2£‘  1 2f  (Q  ')  £ (Q) 

g g g 


(5) 


where 


l-kpJQ')-,173 

£g(°'>  = L— % — ] 


and 


(6) 


rkpe(Qh 

yQ>  -MH  • 


1/3 


Therefore, 


vQ(x)  -e  1 4 Jn  , ! — 


= e')X,P 


> 


2£g(Q-)£g(Q)  “i  t. 


But  from  (la) 


. tt  ® -jx  t 

"J  4 _ ~ e P _ 
e ^ — = vn(x) 


p=l  *p 


Consequently, 


vQ(x)  ” / 

Similarly, 


kD 


0W  ” I— ; v 

0 V 2x  f (Q;)  f (Q)  0 

o o 


vn(x)  = vQ(x) 
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(7) 


(8) 


(9) 


(10) 


- 


r (X)  -*  ' V.(X)  = V (X) 

1 L2xf  (Q  ) f (Q)J  1 1 

g g 


u0(x)  - [— ; ] u (x)  = u0(x) 

0 LZxf  (Q  ) f (Q)  ° u 

g g 

b.  Example: 

Generalization  of  the  Full  Formula  to  a Conical  Surface 
The  formula  (II- 101)  generalizes  for  a conical  surface  to: 

2 1 -jkl 

j mcd  ^[sin  9 # sin0  - Id  sin  8 ' sin  0 + ib  cos  6 ' cos  0 J td 

t-2v  , r - / _1e‘jkD 

+ M — ;un(x)-sin9  sinQv^x)! j 

2TT  cos  e ' cos  e “ 0 1 J k(D; 

OBSERVATION 
POINT 


MAGNETIC 

CURRENT 

ELEMENT 


Figure  IV- 1.  Surface  Ray  on  a Cone 

2 2 

In  (13)  on  a conical  surface  the  pattern  factors  sin  0 and  cos  9 have 
been  replaced  by  sinO'  sin9  and  cos  9'  cos  9 (see  Figure  IV-1),  in  order 
to  preserve  reciprocity. 

Samples  of  numerical  results  for  circumferential  slots  on  a conical  surface 
may  be  found  in  [2], 


V..  Conclusions 


This  report  is  a study  of  the  validity  of  Surface  Ray  Methods  for 
application  to  mutual  coupling  in  conformal  arrays.  The  study  consists  of 

an  analytical  part  and  a numerical  part. 

The  classical  GTD  [1]  includes  leading  terms  in  either  polarization 
of'the  high  frequency  asymptotic  expansion.  These  correspond  to  rays  which 
propagate  along  geodesic  trajectories  of  the  unperforated  conducting  array 
surface  with  a near  free  space  velocity,  and  attenuate  exponentially  during 
their  travel.  Furthermore,  because  of  the  ray  spreading  there  is  an  alge- 
braic ray  divergence  coefficient  — . 

JkD 

In  this  report  It  is  shown  that,  in  addition  to  the  surface  ray  terms, 
diffraction  of  surface  rays  must  also  be  included  in  certain  angular  regions. 


These  diffraction  effects  amount  to  the  inclusion  of  terms 


°(..  ~srr) 


(kD) 


in  the 


deep  shadow  or  O 


f 1 


kD) 


u 


in  the  near  zone.  The  significance  of  the  O 


terms  shows  up  clearly  in  the  circumferential  (H-plane)  coupling  which  would 

be  grossly  underestimated  on  the  basis  of  — - — terms.  The  — y terms  in 

vkD  (kD) 

the  near  zone  are  also  necessary  for  the  near  neighbor  coupling. 

A basic  difficulty  is  observed  in  the  paraxial  (E-plane)  region  for 
circumferential  slots.  It  has  to  do  with  the  break-down  of  the  Airy  function 
approximation  of  H^^(z)  for  small  values  of  the  argument  z.  A correct 
treatment  of  the  dyadic  Green's  function  in  this  region  requires  a different 
approach  than  used  in  this  report.  For  example  a harmonic  series  could  be 
used.  Here,  the  Fourier  integral  with  respect  to  the  axial  wavenumber 
would  be  evaluated  asymptotically  to  second  order  of  the  inverse  distance 
from  the  source.  This  representation  which  is  not  a surface  ray  representation 
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should  yield  a smooth,  valid  transition  to  the  surface  ray  representation. 

A further  question  that  arises  is  the  generalization  of  this  transition 
function  to  the  conical  geometry.  , 

Instead  of  proceeding  along  these  lines  this  report  establishes 
first  the  correct  asymptotic  expansion  of  the  surface  for  circumferential 
magnetic  current  element  on  a circular  cylinder.  This  expression  is  valid 
in  the  circumferential  plane  and  its  angular  neighborhood.  An  approximate 
transition  function  is  then  constructed  that  blends  with  the  above  representa- 
tion in  the  circumferential  plane,  is  finite  along  the  axis  (E-plane)  and 

reduces  to  the  planar  Green's  function  to  of -x-  j in  the  proximity  of  the 

ikDj 

source.  A variant  of  this  formula  is  also  provided 

The  correlation  of  the  numerical  results  in  Section  III  with  those 
obtained  on  the  basis  of  harmonic  series  for  cylinders  with  ka  « 9.  5 and 
19.6  justifies  the  use  of  these  approximate  formulae.  These  above  approxi- 
mate formulae  are  easily  generalized  to  conical  surfaces.  > 

As  mentioned^the  numerical  data  for  between  circumferential 

slots  indicate  a good  agreement  with  the  harmonic  series  results  for  cylin- 
ders with  ka  « 9.  5 and  an  improved  correlation  for  ka  » 19.  6. 

The  agreement  for  Y ^ between  axial  slots  on  a cylinder  with 
ka  « 9.  5 and  based  on  formula  (11-90)  is  very  good,  although  no  phase  com- 
parison is  available. 
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Appendix  A.  Review  of  the  Theory  of  Alternative  Representations 

The  procedure  for  obtaining  electromagnetic  field  solutions  in 
cylindrical  regions  of  uniform  cross-section,  utilize  various  represen- 
tations of  fields  and  their  sources  in  terms  of  a complete  set  of  vector 
eigen  functions,  with  the  resultant  scalarization  of  the  relevant  vector 


problem. 


The  scalar  potentials  are  defined  via  [Ref.  6,  Eqs.  24a  and  25a, 


pp.  196-197] 


-(/)(*) 


, r *}#)(')(P)  (P#) 


(£.£  ) = L 


gz(z.  2 ) • 


In  (1)  the  prime  and  double  prime  denote  TM-  and  TE-  polarizations, 
respectively,  _p  and  _£*  are  the  transverse  position  vectors  of  the  observa- 
tion and  of  the  source  point,  and  g (z,z')  is  given  in  (A18).  The  sub- 
script  i stands  for  a double  index,  one  corresponding  to  the  radial,  the 
other  to  the  angular  domain.  The  mode  potentials  ty/(_£)  and  ff(_p)  indi- 
vidually satisfy  in  the  cross-section  the  Helmholtz  equation 

^/)(/,)(P,w)  + k</)(//)2  t{')(')(p,®)  = 0 ( 


with  periodicity  of  2tt  in  c o of  the  function  and  its  cp  derivative,  and  with 

♦/  (a,  ®)  = 0,  To  ( P,  «P)  I =0 

p=a 

In  (1)  and  (2),  k.  denote  the  radial  eigenvalues  (wavenumbers), 
i 

The  normalization  is  such  that 
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(4) 


TP  |e[/)(")(P)|2  dS=  1, 

cros  s-section 

7 */(p)  * 

where  e. '(P)  = - and  e"(p)  = xz  (5) 

—1  — « / 1 V ' 

\ 

To  facilitate  conversion  of  the  various  alternative  representations, 
it  is  useful  to  define  three  one-dimensional  characteristic  Green's  functions 
as  follows: 

a.  Angular  Green's  Function 


are  shown  in  Figure  A-l. 


Figure  A-l.  Complex  X^-plane  Figure  A- 2.  Complex  X^-plane 


a = arg  (Xu).  3 = arg  (k2  - X^) 


8p(0’p'Xu;Xv): 


pJ — — 

T^V'  U 0W(2)f 


QH(^.'k2-X  a ^ 

VX  "■  U 

^ v 


H(lLf,/k^-X  p ) 

VX  '»  u >y 


where  Q = 1 for  a TM-  case,  Q = for  a TE-  case,  and  p>  denote  the 

< 

greater  and  the  lesser  of  P and  p7  , respectively. 

There  is  no  branch  point  at  X =0,  since,  in  view  of  the  relation 

H(_a2)(Z)=  ei  janH^ ’2\Z)  , 

g_  is  an  even  function  of  J X . 

6 P v v 

It  is  convenient  to  introduce  the  transformations  = v and 
,/k2  - X^  = k^,  whereupon  gp  becomes 

. J Mv  QH(1,(k.a)  m 1 m 

gp(P.  o';Xu;Xv)  = -f  |HW(kt  0<)  - H'2,(kt  P^h'jX  P>)  . 


The  singularities  of  g^  in  the  X^-plane  are  simple  poles  at  the  zeros 

(2) 

of  QHV  (k  a).  They  are  indicated  in  Figure  A-l,  which  also  shows  the 

v v 

contour  of  integration  enclosing  these  singvilarities  in  the  completeness 
relation  [Ref.  6,  Eq.  3-4-98,  p.  326] 

QHv){kta) 

“2Sjf  dXv«p(0'  p/;kt2;Xv)=-jI  I QH^Vka) 

cv  u\  ^yt  V t 


= - p'6  (p.p') 


( 


r 





In  the  X^-plane  gQ  exhibits  a branch  point  at  X^  = k . The  choice 
“2 

of  the  proper  branch  of  ,/k  - X^  is  such  as  to  satisfy  radiation  condition 

for  P “*  ®.  A convenient  branch  cut  and  the  path  of  integration  Cp  for  the 
completeness  relation  [Ref.  6,  Eq.  99,  p.  327] 


rij*  dXugo(p,p':Xu;v2)  = ii^-) 


(17) 


,U 


is  shown  in  Figure  A-2. 

c.  The  Axial  Green's  Function 
This  Green' s function  satisfies 

0^2  + Xu)gz  {2’Z';V  = - 6(z-z') 

subject  to  radiation  condition  for  |s|  — ®.  Therefore 


a (z,  z ) = — 

1 “ 


z-z 


1x11  — °* 


(18) 


(19) 


One  observes  that  gp  has  a branch  point  singularity  at  X^  = 0.  The 
proper  sheet  of  ./X~  is  defined  as  0 > arg  (X^)  > - 2tt  . For  convenience,  if 
the  positive  real  X^-axis  is  chosen  as  the  branch  cut  (see  Figure  A-2),  one 
obtains  the  completeness  relation  for  the  z-domain 


d Xy  gz(z,  z';Xu)  = - 5(z-z')  (20) 

C 

z 

To  facilitate  the  conversion  between  alternative  representations  one 
transforms  the  modal  potentials  Sx  and  S*  of  (Al)  into  the  integral  representa- 
tion of  (21a)  involving  the  above  three  one  dimensional  Green's  functions. 

This  representation,^  turn, enables  one  to  deduce  the  various  alternative 
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representations  merely  by  deforming  the  paths  of  integration  C or  C 

05  Z 

into  Cp;  or  C^.  respectively.  Examples  of  alternative  representations 


for  S are: 


} dXu^r  8o(c-  p'!Xu;XvV®-’’'=xv,«»1*-  *'*v 

c 


(21a) 


u 


^L'L)  = \ y-  2,  cvcos  r <*X  -y g (P,  o';\  ;v2)g  (z,  z';X  ) (21b) 

1 2TT  j v=0  C k -X-  2 


u 


cos 


V,  (IT-  loo-CD7  I) 


1=1 


sin  v^tt 


dXu  ^3-H!ft’(A2-xu 


QH 


(2 


^/k^T  a; 


gz(z,  2';Xu) 


(21c) 


where  Q=1  for  S=S  ' and  Q = for  S=S"  . 

Representation  of  the  potentials  S'  and  S"  in  (Al)  follows  from  the 
relation  [Ref.  6,  p.  198  Eqs.  (32b),  (33b)] 

- V2  S(/)(#)(r,r/)  = G(/)(')(r,r/)  , 


(22) 


where 


G(,)(#)(r,r/)  =^i(/)(")(P)  (f|/,{")*(i/)gz(z,z/)  . 


(23) 


In  a similar  fashion  (21a)  follows  from  [Ref.  6,  p.  285  Eq.  37]  and  (A-22). 


To  obtain  (21b),  the  integral  in  (21a),  with  contour  C ,is  evaluated  as 

05 


a sum  of  residues  at  the  simple  poles  of  g^.  To  derive  (21c),  the  contour 


^CD  *S  ^e^orme<^  *nto  Cp  an<i  the  resulting  integral  is  expressed  as  a residue 
series  at  the  simple  poles  v l o£  8p- 


1 
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Appendix  B:  Curvature  Terms  and  their  Derivatives 


Basically,  there  are  two  curvature  terms,  hard  and  soft. 

r -M  Hv2)(z) 

They  are  derived  from  the  integrals  e 1 — - dv  and 


H(2)  (Z) 

v 

r -jv*  J 

! e — t-t dv  . 

*V(Z) 

a.  The  hard  curvature  term  and  its  derivatives 

For  large  and  almost  equal  argument  and  order  the  asymptotic 

(2) 

expansion  of  the  Hankel  function  (Z)  and  its  derivative  with  respect  to 
argument  is  found  in  the  form  of  a combination  of  the  Airy  functions[4].  The 
leading  term  of  this  expansion  is 


h[2)(Z)  - -J-f-1  w (t) 
H^2)(Z) Lf-2  w'(t) 


for  v » Z and  Z » 


where 


_ v-z  ( zy 

- f - KtJ 


Here 


w2(t)  = Jn  [B. (t)  - j A.(t)]  , 

where  A.(t)  and  B^(t)  are  the  two  independent  real  solutions  of  the  Airy 
differential  equation.  Therefore  w^ft)  satisfies 


2 

( ~T  ‘ 0 w2(t)  = °* 


The  lowest  order  (in  Z)  asymptotic  term  of  the  ratio  of  Hankel 
functions  is,  from(Bl) 
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1 


H®(Z) 

h«v2>'(Z) 


- - £ 


w2(t) 

w2(t) 


(5) 


If  Z » 1,  this  ratio  exhibits  a peak  value  near  the  poles  of  the 
integrand.  Therefore  the  major  contribution  of  the  Fourier  transform  of 
this  ratio  is  derived  from  the  region  of  the  v-plane  where  v » Z.  For 
Z » 1,  one  may  use  the  asymptotic  form  of  this  ratio.  Thus, 


f i: -M  H*>(z) 
-«  H^2)  (Z) 


dv  - -e‘jZ}f2  ’e"jfft  dt 

T w2(t) 


(6) 


where  * > 0,  is  a parameter,  and  the  integration  contour  T depends  on  Z. 

w (t) 

one  notes  that  for  |t  |**®,  is  algebraic  [Ref.  6,  pp.  394-395].  When 

w'(t) 

arg  Z=0,  argt=0  and  t runs  from  -®  to  +°°  along  the  real  axis  in 


the  t-plane,  the  valley  regions  being  the  entire  lower  half  t-plane.  But 

, „ tt  . 5n 

when  arg  Z = - ^ , arg  t = - -g-  a 

with  the  valley  regions  below  the 


, _ TT  ,.5TT  . TT 

when  arg  Z = - ^ , arg  t = - -g-  as  v - ®,  and  arg  t = g-  as  v -*  + <*> 


Figure  B-l.  Complex  t-plane  and  the  Integration 
Contours  for  arg  Z = 0 and  - . 
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contour  T.  Figure  B-l  shows  the  contours  t for  both  arg  Z=0  and  - . 

To  find  a contour  valid  for  either  case,  one  deforms  that  part  of 

the  straight  line  with  Imt  > 0 (for  arg  Z=  -j)  to  the  real  t axis  and,  in 

this  fashion,  recovers  the  usual, broken  contour  T of  Figure  II— 3 running 
. 5n 

f -JT  „ 

from  ® e to  +•  . 

The  hard  curvature  term  Vq(x)  is  defined  by 


VQ(x) 


= re-jxt^dt 

T w2(t) 


where 


X * f * . 


There  are  two  expressions  for  Vq(x);  one  is  for  large,  the  other 
for  small  values  of  x.  For  large  values  of  x,  one  may  close  the  contour  T 
in  the  respective  valley  regions  (i.e.  those  of  arg  Z= 0 and  - at  infinity  in 
the  lower  half  plane.  This  procedure  results  in  a residue  series  evaluated  at 
the  simple  poles  of  w^(t)  = 0: 


jx  t w_(t  ) 


Z-JX.  t W - 

e P“T 

P=1  W2 


~2<V 


where  t = |t  | e J 3,with  Ai  (-  1 1 |)  = 0.  Numerical  values  of  |t  | may 


P P 
be  found  in  [7], 


From  (B4),  one  finds 


* -jxt 


V*>  = - I V 


p=l  p 


For  large  values  of  x,  this  residue  series  converges  rapidly.  In 
contrast,  for  small  values  of  x the  residue  series  becomes  poorly  convergent. 
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An  alternative,  more  efficient  representation  in  this  range  is  a "power 

series"  in  x which  can  be  established  by  an  asymptotic  expansion  for 

w2(t) 

large  values  of  t of  and  a term  by  term  integration.  To  this  end, 

w2(t) 

one  may  rotate  the  contour  T clockwise  by  90°.  For  arg  Z=0,  T maps 
into  a Bromwich  contour  t ' by  letting 


t = j s. 


whereupon 


„ w (is) 

vow =j  r e*s-r-r ds 

-n  ws(js| 

od  e 


For  arg  Z - - — one  rotates  the  contour  clockwise  by  120°  by  the 
transformation 

. TT 

. . 

t = j e s , 

which  again  yields  the  same  Bromwich  contour  T i.  e. 

. tt 

.TT  , /.  n \ 

x is  w2 v e s;J 

V (x)  = j e ' “ Aa 


v*o 


Following  [8],  one  has 


W2(js)  'j4T  1 . 7|  . 21 

—TT 7~e  L“f • TT  + ,,7/2  + “7V  5 + ° 

w^ljs)  s'  4s  32s  /-la- 


s' 4 s 32  s 


VI 64 1 


(°-13/2;] 


Similarly,  one  finds 
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where  I"*  denotes  the  inverse  Laplace  transform.  Therefore  for 
arg  Z=0,  one  finds 


VQ(x)  = 2tt  rl 


w2(js) 

w2(js) 


+ JJL-  . . ) 

512./I 


and  for  arg  Z = - ^ , one  obtains 


(17) 


Note  that  (B17)  and  (B18)  have  the  same  form,  but  the  former  is  valid  for 


real  x while  the  latter  for  x complex.  In  either  case,  if  one  defines  Vq(x) 
as  the  "power  series"  shown  in  the  bracket  of  (B17)  or  (B18) 


Defining  v^(x)  as 


, +ipL  £ *3  .JL 


¥ *9/2+. 


/r  64 


(x)  ={ 


-J  Tri  3/2  V"  - JX  t 

2y^  X V J p 


TT 

-Jz 


4 p=  1 


for  small  x 


for  large  x 


(23) 


(24) 

I 


one  has  Vy(x)  = - 2TTj 


. e 


. TT 

-J4 


172  vi(x)  * 


(25) 


7.J*  x 

For  ten  residue  terms  and  four  power  series  terms,  the  difference 
between  (B25)  and  (B26)  at  x = 0.  6 is  about  0.42  percent  in  the  amplitude 
and  0.  126  degree  in  the  phase.  Figure  II- 5 shows  the  cross-over  for  v^(x). 
The  second  derivative  of  the  hard  curvature  term  is 
. w?(t) 

V'(*)  = - IW3*'-?— dt=  - V <*)  , 

T w2W 


(26) 


With 


•(x)=< 


one  has 


i + 41  x3  + *9/2  +.... 

12  Jj  52 


4 M H 5/2 y — -ix,P 

3^e  x £tpe 

D = 1 


. tt 

“JZ 

v2(*>  - - hfpr-  v2<*> 


for  small  x 


for  large  x 


(27) 


(28) 


(29) 


101 


The  cross-over  point  is  in  the  proximity  of  x = 0.  6 where  the 
difference  between  (B29)  and  (B30)  is  about  1.7  percent  in  the  amplitude 
and  0.  101  degree  in  phase. 

b.  The  Soft  Curvature  Term 

r -iv«  H<v  ’ ,(z) . -JZ»  r -txt  w2(l) 

I*  dv  J'  ^ 

■ 89  \J  ' T 

where  x = f$,  and  t is  the  same  contour  as  shown  in  Figure  B-l. 

The  soft  curvature  term  Uq(x)  is  defined  by 

r -ixtw2(t) 

°ow  * J • 1 soiy dt 

T 


In  an  analogous  fashion,  the  residue  series  representation  of 


UQ{x)  is 


Z-jxt 
e p 


where  t = It  | e 3 , Ai  (-  |t  |)  = 0.  This  series  (B34)  converges  rapidly 
P P P 

for  large  x. 

Similarly  for  small  x,  with  t = js,  one  has  (since  arg  Z * - * will 
not  generate  a new  function,  we  consider  only  the  case  arg  Z=0  here) 

wj(Js)  Ls?  4s2  32s7/Z  64  s2  V 

Hence. 


i. 


Appendix  C.  Asymptotic  Series  for  Laplace  Type  Integrals 

Derivation  of  the  complete  asymptotic  series  for  the  case  of  an 
isolated  saddle  point  for  the  following  canonical  integral  is  reviewed  below. 

I(fi)  = f enq(a)  F(a)  da  (1) 

SDP 

where  » 1,  SDP  denotes  the  steepest  descent  path  while  F and  q are 
analytic  functions  of  a. 

If  q(a)  has  a first  order  saddle  point,  i.  e.  q 7(a  ) = 0 but 

s 

q"(as)  4 0.  then  one  transforms  q(tt)  into  a simpler  function  T(s)  via 

q(«)  = T(s)  = - j - s2  (2) 


With  this  transformation  the  saddle  point  maps  into  the  origin  of  the 
s -plane.  Therefore 

q(*s>  = - j • 

Consequently,  the  integral  becomes 

nq«V  r .ns2 

I(fi)  = e 3 j G(s)  e’us  ds  , 


(3) 


(4) 


where 

G(s)  = F(a)  ||  . (5) 

The  proper  branch  of  s = ,/-j  -q(a)  is  determined  by  the  orientation 
of  the  SDP,  i.e.  by  the  phase  angle  of  ||  at  the  saddle  point.  One  then 
expands  G(s)  into  a power  series  about  s = 0,  and  evaluates  the  integral  in 
'4)  by  carrying  out  the  integration  term  by  term.  The  radius  of  convergence 
of  this  power  series  will  be  determined  by  the  numerical  distance,  or  the 
distance  from  s = 0 to  the  nearest  singularity  of  G(s).  The  numerical  distance 
will  determine  the  error  of  the  asymptotic  expansion.  Thus, 
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C*  OD 

Kn)-e  8 J LG(0)  + sG'(°)  +|rG'(°)  +•  ..]e‘fis  ds 


c(2n)(o) 

" L _/2  n\t 


T2iTl  0 Yn+1 

n=0  ““T" 


(6) 


where  ~(Z)  denotes  the  Gamma  function  and  G^2n^(s)  the  2nth  derivative  of 
G(s).  Note,  odd  terms  do  not  appear  since 


l'  n -CJs  ds  = 0 for  n odd. 
s 6 


(7) 


Figure  D-l.  Surface  Ray  on  a Ground  Plane 


»<£.£•  = - iVn-L6'  6 (>  - i IE"  ^f) 


-jkD 


(kD) 


kD 


and 


/4  . k2  w / « (,  . I \ e'jkD 

Bl»£>  = - JlrrM*S  n \ “ ^"icD-  J kD~" 

Here  M denotes  an  arbitrarily  oriented  magnetic  point  source  and  r has 
been  replaced  by  D the  distance  between  the  source  and  the  observation 
point. 


Formulae  (8)  and  (9)  are  in  a vector  form  and  are  thus  invariant  to 


the  choice  of  coordinate  system. 


Appendix  E.  and  for  Circumferential  Slots 


This  program  calculates  the  mutual  admittance  Y12  (and  S^) 


between  "single  mode"  circumferential  rectangular  slots  on  a Conducting 


Circular  Cylinder.  The  slot  illumination  is  TE.n. 


The  geometry  is  simply  illustrated  on  a developed  cylinder. 

Z 


ZSEP  = (z-z') 


nd 

HE  2— SLOTS 


a<p-QZ  (IPi) 


REFERENCE  SLOT' 


Figure  E-l.  Developed  Cylinder-Circumferential  Slots. 


In  Figure  E-l  the  center  of  the  reference  slot  is  fixed  at  r = a, 
CO / = 0,  z ' - 0 of  a standard  cylindrical  coordinate  system.  Thus,  the 


reference  slot  is  located  at 


Q1  (1)  = Radius  of  the  cylinder, 


Q1  (2)  = cp#  = 0°, 


Q1  (3)  = z = 0. 


The  locations  of  the  centers  of  the  other  slots  are  generated  by  DO 


LOOP  56.  Each  step  of  increase  in  is  10°  for  R=  1.991"  and  10.133° 


for  R = 3.777",  For  different  steps  value  the  step  size  in  DO  LOOP  56 


should  be  modified.  The  symbols  are 


Q2  (I)  * radius  of  the  cylinder 


Q2  (Ipl)  * *o  • «c„oter 
02  (Ip2>  , 


The  long  side  of  the  slot  is  divided  into  N subdivisions  and  the 
short  side  M subdivisions.  Thus,  the  total  number  of  sub-areas  in  one 

t 

aperture  is  M x N. 

/ c \ 

The  Green's  function  of  the  surface  magnetic  field  \^H j is  found 
in  the  subroutine  "Fock",  where 

HTD  = . 

If  ITT  = 1,  denotes  (or  S^)  calculated  by  the  "variant1'  formula. 

If  ITT  4 1,  the  "Full  Formula"  is  used. 

In  the  program: 

Yjj  is  the  self-admittance  of  the  aperture  in  TE^q  illumination, 

Yg  - the  characteristic  admittance  of  the  feeding  waveguide, 

ZSEP  - the  axial  separation  of  the  reference  slot  and  the  second 

slot, 

WL  - the  wavelength. 

The  curvature  terms  Vq(x),  v^(x)  and  Uq(x)  are  termed  in  the  program 
VV,  VV1  and  UU,  respectively.  X(=  f^f)  is  called  XX  in  subroutine  "Fock". 
Aperture  size  is  a x b . 

The  program  is  composed  by  three  parts;  (1)  the  main  part, 

(2)  the  subroutine  SLOCA  (which  treats  the  locations  of  the  subdivisions  of 
the  second  slot.)  and  (3)  the  subroutine  Fock. 
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V G LEVEL  21  MAIN  HATE  = 76190  13/14/21 

C MUTUAL  ADMITTANCE  OF  Cl  RCUM  F ER  CT_IA  L SLOTS  ON  A CYLINDER 

COMPLEX  HT,HR,YY (200) ,SM,Y1 ,CEXP 
01  Mi  NS  I ON  02(200), 01(3) 

REAL  K,LO 

COMMON  /AREA1/ITT,N,M,PI  ,AA, Bb,CC, K , WL , 0 EL , SX , S Y 
C Ql(i)  = *,  Q1(2)*PMI,  Ql(  3 ) = Z. 

C IN  THE  INPUT  WL  IS  IN  CMS,  ALL  OTHERS  ARE  IN  INCHES. 

RE  AD  (5,  <,11I01<  l)  ,Q1<2)  ,01(3)  ,WL,ZSEP 
READ (5,412 )YG»  Y l 

411  FORMAT ( 5FI C.5 ) 

412  FORMAT ( 3E 1 5. 5) 

ITT*  l 

C IF  I TT= l DENOTES  THE  FORMULA  OF  NO  VI  AND  DROPPING  VO  FROM  XC.3. 

ki=Ol(  l)*2.54 
Pl=3. 14159 
TP=2. 33811 
TFB* l. 018  79 
NDl*  19*6 
ND2  *ND  1*2 
ND 3* NO  1*3 

C Yl=S  ELF  ADMITTANCE,  YG*ChARACTFR I S T I C ADMITTANCE 

Y1=Y1*YG 

YA=SQRT(REAL(Y1)**2+AIMAG(Y1)**2) 

YP  = ATAN2( AIMAG(Yl) , REAL  ( Y 1 ) )*130./P! 

C WAVEGUIDE  SlZE=AA"vB8",  WL  = WAVEL  ENGTH  IN  CM. 

AA* .9*2.54 
BB  = .4*2.54 
K=2.*PI/WL 
SSX* AA 
N=  14 
M=2 

l\NM=N*M 

WP. I TE(  6, 35  ) NNM 

35  F0RMAT(1HI,1X,I3,1X, 'POINTS  IN  THE  APERTURE,  CIRCUM.  SLOTS') 

C GENERATE  THE  ARRAY  OF  ThE  SECOND  SLOTS. 

WR 1 T E(  6 » 22 14  ) 

2214  FORMAT ( 20X , • K.  CYL.'»l6X,'Z  IN  CM.',  3X,'PHI.  IN  DEG.') 

DO  56  1=1, NDl, 3 
IPi*  1+1 
IP2=I*2 
11*1  P2/3 
02(1  ) = 01(1) 

Q2(IP1 >=10.*( IP2/3-1 )*P 1/160. 

02(IP2)=ZSEP*2.54 

IF( I. GT. NO  1/2)  Q2( I P 1 ) =-2.*P I + 02 ( I Pi-NDl/2) 

0P=Q2 (IPl)*180./PI 

56  WR 1 T E( 6 , 82 ) 1 1 , 02 ( I ) ,02 ( I FI ), C2 ( I P2 ), QP 

KRITE(6,81)  (Q1(L),LS1,3) 

WRIT  E ( 6 ,93 ) N,M,AA,B6,WL 

81  FORMAT (IX, 'R1=',F12.5,10X, 'PH  l* • , F12.5 , 10X, • TH  1*',F12.5) 

82  FORMAT (IX,  !3«10X,4F12«5) 

83  FORMAT ( IX, 'N*' , 1 3 , 10X , • M =• , I 3 , /I X , • A* • , F 12.5 , • CMS ' , 10X , 
l'B*' ,F 12.5, 'CMS' ,10X, 'WAVELENGTH*' ,F 12. 5 , ' CMS ' ) 
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V ft  LEVEL  21 


MAIN 


DATE  = 7 6 190 


13/14/2  L 


WRIT  El o » 97  6 ) Y1,YA,YP,YG _ 

976  FORMAT (IX,  • SCLF-AOMI T= • , 2E2C  .6  ♦ l OX  , • MAC,*  • , F 1 5. 5 , 10X,  ' PH= ' , F L 5. 5 , / 
11X,'CHAC.  ADMIT.** ,E20.5) 

P01=01(2)*PI/180. 

SX=A A/ N 
SY*B8/M 
DEL  = SX*SY 
WRITS(6,l9l) 

191  FORMAT ( 1H1 , IX. *Y  DENOTES  MUTUAL  ADMITTANCE,  S ISOLATION') 

C NUMERICAL  INTEGRATION  BEGINS 

WRITE! 6,547) 

547  FORMAT  UX,  'PHI  • ,11X,  'REAL'  , 18X,  ' IMA*  ,14X,  'MAG'  ,16X,'  PHASE* , 15X, 

1 * DB ' ) 

DC  50  I2*l,NDl,3 
IJK*(I2-l)/3*l 
YY (I JK )*(!., 1.) 

IF(  I JK.GT.  13.AND.IJK.LT.3DG0  TO  50 

I F ( UK.GT.36  JGO  TO  50 

IP  1 = 12* l 

IP2=  1 2 2 

XQ2-Q2 (IP1)*RI 

Y02=Q2(IP2) 

HT  = ( 0 • » 0 • ) 

HR=(0. ,0.) 

I F ( YQ2. NE. 0. 0 ) GO  TO  1074 
JFf 1P2/3.LT.4) GO  TO  50 
1074  CONTINUE 

DO  41  I =1, N 

snx=  (float (n+1  )/2. -float  in  )*sx 

E2=C0$(PI*SNX/(SX*N) ) 

X(J2C*XU2  + SNX 
DC  41  J*l,M 

SMY=(FL0AT(M+l)/2.-  FLCAT(J) ) ^SY 
YQ2C=Y02-SMY 

41  CALL  SLOCA ( E2 , R1 , XQ2C , YQ2C , HT , HR  ) 

YY  ( I JK  ) =-HT*2./(AA*B6) 

I F ( I JK.GT. 12 ) GO  TO  5C 
IJ  = 10*(I JK-l) 

TYY=SOKT(REAL(YY(IJK  ) )**2+AlMAG(YY(  IJK  ))**2) 

TYP*ATAN2( AIMAG( YY( IJK  ) ) , RE AL ( YY ( I JK  )))*180./PI 
TYL=20.*ALCG10(TYY/  YA) 

WR1T£(6,958) IJ  , YY(IJK  ) , T YY , Typ , T YL 
SM=-2.*YG*YY(  I JK)/(  ( YG*-Y1  + YY<  IJK  ) ) *( YG*Y l- YY ( IJK)  ) ) 
SA=S0RT(REAL(SM)**2+AIMAG(SM)**2) 

SP=ATAN2  ( A IriAG  ( SM)  , REAL  ( SM)  )*18U./PI 
SL=20.*ALGG10(SA) 

WRITc(6,959)SM,SA,SP,SL 
50  CONTINUE 
NN*N01/o 
DG  31  I J * l ,NN 
IFC  I J.LT.13)G0  TO  31 
I JK* 10* ( I J-l ) 

YY ( I J)*YY( I J ) +YY ( I J+NN  ) 
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SM=-2.  *YG*  YY 1 1 J J / L(  YG«-  Y 1 +Y  Y ( I J II  * ( YG  ♦ Yl-YY  U J 1 1 1 

SA»SU«T ( REALISM )* *2 >AIM ACC SM)**2  I 
SP=A  TAN2  ( A I MAG  ( SM)  ,KEAL  ( SMJJ  *180./  <>I 
SL«20.*ALOGlO( SAI 

YAM  SORT  (AIMAGCYYC  I J I ) * *2*RE  AL  C Y Y ( I J ) )**2) 

YA  WM*  YAM/YG 

YAR=AT AN2C  AIMAGC  YY ( I J ) ) , RE AL < Y Y ( I J ))  ) * 1 8 0. /P 1 
YA  RR  *-360. +YAR 
TYL»20.*AL0G10( YAM/  YA) 

TYPD*YAR-  YP 
01  SK  *!• 0 
TYPDD=TYPD-DIS* 

TYOOC=YAR-DISK 

WR1TE(6*958) I JK , Y Y ( I J) ,YAM,YAR,TYL 
WRIT  EC  6*959) SM,SA,SP*SL 
31  CONTINUE 

63  FORMAT ( IX, *S=' ,2E14.5,5X,‘ AMP=» , E14.5.5X. 'PHA=* , EL4. 5,5X, 
1«0K-S=' tEi^.5) 

33  7 FCRMAT  ( 3 X , • YG ' , I 3 , • * • , 2 E20  .6 , 1 0 X , • MAG=  ' , EL5.5,l0X,,oH=',F15.5, 
13X,F  7.  2) 

827  FORMAT (IX, • 08- Y= • , E20. 6 , 10X , • FH-CI F= • ,3F  15. 5 ) 

33  FORMAT  (IX,  *Y0'  ,12 , • *• ,2E20.6  ,/ ) 

958  FORMAT! IX, I3.5X,* Y=* ,2S  1 5.  5 , 5X  , E 1 5 .5 ,5X , FI 5. 5 , 5X , E 15 .5  ) 

959  FORM  AT ( 9X , • S= • , 2E 1 5.5 , 5X, E 1 5. 5 , 5X ,F 15, 5 , 5X , 5 1 5. 5 , / ) 

STOP 

END 
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SUBROUTINE  SLUCA ( E2 » R 1 . X C2C . YU/C « KT  * HP ) 

COMPLEX  HHT,HTD,HTS,HT,HR 
K^AL  KtLU 

COMMON  /AREA1/ITT  , N.M.P!  ,AA,  BtJ.CC.K,  WL,OEL,SXtSV 

DEL=SX*SY 

DO  l l I =1,  N 

SNX=  (FLOAT  (N*-I  )/2. -FLOAT  (I  ) ) *SX 
EC  = C OS (PI*SNX/ (SX*N ) ) 

XO  = XvJ2C— SNX 
DO  ll  J=l,M 

SMY=(FLOAT(M*l)/2.-  FLOAT ( J 1 )*SY 
C SMY=  C J — 2 ) * SY 

Y0=YQ2C«-SMY 
IF (XC.EQ. O.OJDO  TO  21 
TCL=ATAN2( YO.XD) 

GO  TQ  23 
21  TCL*PI/2. 

C IF(YO.LT.O.)TCL=-PI/2. 

23  LQ=K 1/ (COS (TCL 1**2 ) 

FF=EXP ( ALUGt  K*Kl*AflS(COS  (TCI  ) ) / 2 » ) /3  • ), 

ANG=XO/Rl 
XX=FF*A3S ( ANG) 

SOK=SQkT { XD**2  + YU**2  ) 

CALL  FOCKl  ITT,K,Rl,CCtXX,Pl  , FF  , LO  , SOK , TC  L , HT  0 , HTS  ) 
HHT-HTD+HTS 

HT=HMT,"tC*F2*DEL**2  + HT 
HR= (0,,J,  ) 

11  CONTINUE 
RETURN 
END 
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SUBR OUTlNE-.EOCKt  III,  K ,R1  ,CC.,  XX , PJ.,  £Ej  LU  ,.DIS_,I ,ilLD.,-HI£L) 

COMPLEX  H4,H5 

COMPLEX  VV ,LU,VV0,VP ,CEXP,HTO,HTS, H3,H1, VV1, VV2 
COMPLEX' KKDTP,B3B,AA1 , A A2»TT0, YYZZ ,K KDT  f KKD , TT0P,CQN JG 
REAL  K,LQ 

DIMENSION  TFB(10)»TF(10) 

C TFB  ARE  THE  ROOTS  OF  THE  AIRY  FC T DERIVATIVE,  TF  THE  ROOTS  OF  AIRY  FCT, 

DATA  T FB/ 1.0 13 19, 3. 24819,4.820 10, 6. 16331,7. 37218,8.43849, 

19. 53545, 10.5277, 11.4751 ,12.3848/ 

OATA  TF/ 2. 33 8 11, 4. 08 795, 5. 52056,6. 78671,7.94413, 

19.02265, 10.0402, 11.0085, 11 .9 360, 12 .8288/ 

PI=3. 14159 
PRO=SQRT (PI )/2. 

ALG=RL/COS (T ) 

TLG=R1/$IN (T)**2 
H5=C  EXP ( (0.,1. )*PI/4.) 

H8  = SURT  (f>I  ) 

Y=l./377. 

YYY=K*DIS/ (2.*£XF( AL0G(K«Rl/2. »*2./3.) > 

IF  ( XX. LT. 0.8000  ) GO  TC  53 

VVD=<0. , 1.  )*PI /3. 

VP=(0.,-i. )*XX*CEXP( (0.,-l.)*Pl/3.  » 

VV=(0.,0.) 

VVl=VV 
VV2  = VV 
uu=vv 

RESIDUE  SERIES 
00  3 7 1 = 1,  10 

VV=CEXP«  VP*TFB(  I ) + WD)/TFBl  ! ) +VV 

VV1=CEXP(VP*TPB( I)  J+VVl 

VV2  = TFB(  I )'  *C  EXP  { V P*T  F B ( I )-VV0H-VV2 
37  UU=CEXP(  VP*TF(I)  J-MJU 

VV=VV*SGRT(P1=XX)*CEXP( ( C. ,- l. > *P I /4. ) 

VVl=VVl*2.*S0RT  <P1*XX**3  >*CE XP ( (0 . , 1 . ) *P 1/4. ) 

VV2= VV2*4. *(0. ,1.  )*CEXP( ( 0. , 1 . ) *P I /4. ) *S QRT ( P I*XX**5 )/3 . 
UU=UU*2.*SQRT(PI*XX**3)*CEXP( (0.  , i.)*PI/4.  ) 

GO  TO  54 
POWER  SERIES 

53  VV= ( l.-CEX  P( (0. , 1. )*P I/4.)*SQRT( PI *XX**3  J/4. 
l + 7.*(0.,l.  )*XX**3/60.+7.*S0RT(PI*XX**9)*CEXP(  (0.,-l.  )*PI/4»  )/5l2.  ) 

VVl= l.+CEXP{ (0., 1. )*P 1/4.)* SORT (PI *XX«*3 ) S2.-7.* (0. , 1. ) 
1*XX**3/12.~7.*CEXP( ( 0 • , - 1. ) *P 1/4 . ) *SQRT ( PI *XX**9 ) / 64 . 

VV2=l.+7.*(0.,  1.  )*XX**3/12.+7.*CFXP(  (0.,-l.  )*PI/4.  )*SQRT  ( PI  *XX**9 ) 

1/32. 

UU=  ( l.-CEX  P(  (0.,  1.  )*PI/4.)*SCRT (PI *XX**3)/2. 

1 + 5.* (0. , 1. )*XX**3/12.+CEXP( (0. ,-l.  )*Pt/4. )*5.* 

1SQRT (PI*XX**9)/64.  ) 

H4  = l .-(l.+  SINU  )**2  )*H5*SQRT  CP  I*YYY**3)/2.+(  5.  + 7.*STN(  T )**2) 

1*(0. , l. )*YYY**3*C0S(T)**2/l2.+(5.+7.*SlN(T)**2)*SQRT(P!*YYY 
1**9)*CQS(T)**4/(64.*H5) 

54  CONTINUE 

H1=(K**2*Y/(2.*PI ) )*CE  XP ( ( 0. ,-l. ) * K*  01  $ ) / ( K*D I S ) 

HTD=(0.,-1.)*VV*(SIN(T  >**2-2*SIN(T  j* *2*( 0. , 1 . ) / ( K*OI S ) +COS( T ) **2 
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‘ C MUTUAL  COUPLING  OF  AXIAL  SLOTS  ON A CYLINDER 

COMPLEX  HT  ,HR  ,YY(200),SM,Yl,CMPLX 
DIMENSION  C2I200) »Ql (3) 

REAL  K ,Ls) _ . 

COMMCN  /AREA1/N, Mv PI , AA, EB ,CC, K , WL , DEL ,SX , SY 
C IN  THE  INPUT  WL  IS  IN  CMS,  ALL  OTHERS  ARE  IN  INCHES. 

C Ql(ll-Rt  Q 1 ( 2 ) =PH I , Q1 (3  )=Z. 

READ (5,411 )Ql(l),Qi (2 ) , QK3  ) ,WL , ZSEP 
REA0(5,412)YG. Yl 

411  FORMAT ( 5F1 0.5  > 

412  FORMAT ( 3E1 5. 5) 

Rl=Ql(l)*2.54 

P1=3.14I59 

ND 1=19*6 
N02  =NC1*2 

ND3=NO 1*3 ‘ 

Y l=Y  1*YG 

YA=SQRT(REAL!Y1)**2*AIMAG(YD**2) 

YP  = ATAN2(  AI MAGI  Yl) , REAL  (YU  1*180. /PI 
C'  APERTURE  STZE=AA"*8B" 

AA  =•  4*2.  54 

* BB=.9*2. 54 

K=2.*PI/WL 

SSX=AA 

__ N=2 

M*10 

NNM=N*M 

WRITE!  6, 35)  NNM  

35  "FCRMATUHl  , IX, 13, IX,  •PCIATS  IN  THE  APERTURE,  AXIAL  SLOTS') 

C GENERATE  ThE  ARRAY  OF  THE  SECONO  SLOTS. 

WRITEI6, 22141  

2214  FORMAT (2 OX , ' R.  CYL. ' , 16X  2 IN  CM.*,  3X,'PHI.  IN  DEG.') 

DO  56  1=1, NO l, 3 

1PI-U-1 

IP2»1*2 
1 1 * I P2  / 3 
Q2II  )=Q1  ( 1 I 

Q2 1 1 PI ) > 10.* 1 1 P2/3- 1 1 *P I /1 80 . 

Q2( I P2 ) = ZS  EP*2 .54 

IF (I .GT.NOl/2)  C2(IPl)»-2.*PU-Q2(IPl-NDl/2) 

QP*02(IPl) *180./PI 

56  WRITE(6, 82)11, 02 (l)«Q2(IFl)«Q2(IP2),QP 

WRITE! 6, 81  ) I01ILI ,L«l,3)  _ 

WRITE!  6, 83)  N,H,AA,BB,WL 
C 

81  F0RMAT(1X,'R1=*,F12.5,10X, *PH  1= • , F12.5, 10X, *TH  1«*,F12.5) 

82'  FORMAT!  IX,  13, iOX,4Fl2.5) 

83  FORMAT (IX, *N=* ,I3,10X,*M«* ,13, /IX, *A= * , F 12.5, * CMS • ,1  OX, 

1*8='  ,F 12 .5 , 'CMS' ,10X,' WAVELENGTH*' ,F 12. 5, 'CMS ') 

WRIT  6(6,976)  Yl,YA,YP,YG 

976  FORM AT ( 1 X , • SEL F- ADMIT* • ,2E20.6,10X, 'MAG* • , E15.5, 10X, 'PH=' ,F15.5,/ 

1 IX, 'CHAR .—  AO Ml T.*' , E20.5)  

FQ1=01 (2 )*P 1/ ISO. 
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SX=AA/N  

SY*BB/M 
DEL*  SX*SY 

WR  IT  E(  6 , 19  l ) 

191  FORMAT (1H1  , IX, 'Y  OENCTES  MUTUAL  ADMITTANCE,  S ISOLATION'! 

C NUMERICAL  INTEGRATION  BEGINS 

WRIT  E ( 6, 54 7 ) 

547  FORMAT ( IX , 'PHI ' ,1 1 X,  • REAL*  1 18X  , • I MA*  ,14X  MAG*  , 16X,"' PHASE' , 15X, 
l'OB*  ) 

CO  50  12=1  ,N01  ,3  _ 

I JK* (12—11 /3+1 
YY( I JK 1*1 l.,l. ) 

IF (I JK.LT. 3.AND.ZSEP.EC.0. )GC  TO  50 
IFdJK.GT. 19. ANO.IJK.LT. 31IG0  TO  50 
IP1*  12  + 1 v 

IP2=I2+2  

XQ2*Q2( IP1 >*R1 
YQ2=  C2  1 1 P2  ) 

_HT*  ( 0.  ,0. ) 

HR* ( 0. ,0.  ) 

00  41  I*l,N 

SNX*  (FLOAT  (N  + D/2. -FLOAT  (I  J )*SX 

X02C  =XC2  +S  NX 
00  41  J*1 » M 

SMY* (FLOAT (M+l >/2.-  FLOAT(J))*SY 
E2=C0S(PI*SMY/(SY*M) ) 

Y02C *YQ2-SMY 

C SLUCA  TREATS  THE  SUBDIVISIONS  OF  THE  2ND  APERTURE. 

41  " CALL  S LUCA (E2 , Rl,XQ2C, Y Q2C ,HT,HR I 
YY  ( I JK  } *-FT*2 •/ ( AA*8B) 

IF ( IJK.GT.12JG0  TO  50 

TYY=SORT(REAL(YY(IJK  h **2+A IMAG (YY  f I JK  1 )**2  ) 

TYP=ATAN2 ( A I MAG  ( YY  ( I JK  ) >,  RE  AL  ( YY  < I JK  )))*180./PI 

TYL=20.*ALCG10(TYY/  YA) 

I J=10* ( I JK-1 ) 

WR  IT  E(  6,95  8 ) I J , YY ( I JK  ) ,TYY,TYP ,TYL 
SM=-2. *YG* YY ( I JK )/ ( ( YG+Y l+YY ( I JK > ) *( YG+Yl-YYl IJK) ) ) 

SA*SCRT ( RE  AL ( SM)**2*A I M AG( SM )**2 ) 

SP  = AT  AN2 ( A IMAG( SM) , REAL ( SM ) ) *180./PI 
Sl=2  0. *ALOGl 0 ( SA ) 

WR IT  E( fc, 95  5 ) SM ,SA , SP , SL 
50  CONTINUE 

NN-ND1/6  

"bo  3i  ij-i.NN' 

I F ( IJ.LT.13IG0  TO  31 

I JK=10*( IJ-1)  _ 

YYU  J)  *YYl  IJ  )+YY(i  J*NN  ) 

YAM=SQRT  (AIMAG(YY(  I J ) ) =>*2+RE  AL  (YY( IJ  ) )**2) 

YAR=AT  AN2 ( AI MAG ( YY ( I J) ) , REAL (YY (1 J ) ) )*180./PI 

TYL*20.*ALCG10( YAM/  YA)  ‘ ~ 

SM=-2.*YG*YY( I J I / ( (YG+Y14YY( I J ) ) * ( YG+Yl-YY ( IJ))) 
SA=SQRT(REAL(SM)**2«-AIMAG(SM)**2) 

SP*ATAN2(AIMAG(SM),REAL( SM ) )*180./PI 
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SL*20.*AL0G10< SA) 

*RITE(6t958>  I JK.YYU  J)  » Y AM  »YAR»TYL 
VKIT  £( 6 1 95  9 ) SM«  SA»  SP  tSL 

31  CCNT INUE  

'958  FORMATUX,I3,5Xf ‘Y**, 2E  15.5, 5X,E15. 5, 5X,F15. 5* 5X, 6 15. 5) 
959  FORMAT ( 9X , • S» * , 2E 15.5 .5 X, E15.5, 5X , F 15. 5 ,5X, E15.5, / 1 

STOP 

END 
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SUBROUTINE  SlOCA(E2,Al,  XC2Ct_YQ2C_tHT, HR) 

COMPLEX  HHT»hTCfHTStHT»HR 
REAL  K,LQ 

COMMON  /AREAl/N.M, PI,AA,EBtCCf K,WLtOELtSX,SY  

DEL»SX*SY 
00  11  I =1,  N 

SNX* (FLOAT  (N>i)/2.-FL0ATU  ))*SX 

X0«X02C-SNX 

CO  11  J*l t M 

SHY* ( F L OA T ( M ♦ l ) /2 . - FLOAT! J))*SY 

EC*COS I pI*SMY/ (SY*M) ) 

YD*Y  C2C  + SMY 

IFtYC.EQ.O.O  )GQ  T0_2i 

f CL* AT  AN2 ( YC  * XD ) 

GO  TO  23 

21 TCL*  0.0  _ 

IFIXC-LT.O.»TCL*PI 
23  LQ*R 1/ (COS ( TCL )**Z) 

FF*EXP  ( AL0G(K*R1*ABS(C0S  (T  CL  ) ) / 3 . ) 

ANG*XD/R1 
XX  *FF*  ABS! ANG ) 

SOK*SQRT (X  C**2*YD**2  I 

C FOCK  TREATS  THE  CURVATURE"  TERMS  AND  THE  GREEN • S FUNCTICN. 

CALL  FOCK(K,Rl,CC,XX,PI,FF,LO, SOK , TCL ,HTD#HTS I 

HHT*HT  0+-HT  S 

.HT«HHT*EC*E2*DEL‘**2<-HT 
HR* ( 0. »0« ) 

U CONTINUE 

RETURN 

END 
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SUBROUTINE  FOCK( K,R1 ,CC,XX ,PI ,FF,LQ, DI S, T , HTD» HTS) 

COMPLEX  VV  »UU » VVO, VP ,CEX  F»  FTD ,HTS  »H3  ,H1»VV1»VV2*H4*H5 
REAL  K 

DIMENSION  TF6 ( 10) « TF ( 10 ) _ 

C TF B ARE  THE  ROOTS  OF  TFE  AIRY  FCT  DERIVATIVE',  TF  THE  ROOTS  OF  AIRY  FCT. 

CAT A TFB/1 .01879,3.24819,4.82010,6.16331,7.37218,8.48849, 

19. 53545,10.5277,11.4751,12.3848/ 

CATA  TF/ 2. 33811 , 4. 08 795 , 5. 52056,6. 78671,7.94413 , 

19.02265, 10.0402,11.0085,11.9360, 12.8288/ 

PI-3.14159 

F5-CEXPI (0.,l.)*PI/4.) 

H8=SQRT (PI  ) 

Y=l./377. 

YYY=K+DIS/ I2.*EXP(AL0G(K*R 1/2. )*2./3.) ) 

IF  ( XX. LT.  C.8000  ) GO  .TC  53 

VVD= ( 0. , 1 . )*PI/3. 

VP  = <0.,-1.  1*XX*CEXP( (0.,-i.)*PI/3. ) 

VV=IO. ,0.) 

VV  1=  VV 

VV2=VV 

LU=VV 

C RE  SI CUE  SER I ES . 

CO'  3 7 1 = 1,10 

VV=CEXP(VP*TFB( I )+VVC)/TFB( I l+VV 
VV1=CEXPCVP*TFB(I))+VV1 
37  VV2=TF B( I ) *C EXP ( VP*T  FBI  I )— VVO ) +VV2 

VV=VV*SQRT (PI*XX)*CEXP( (0.,-l.)lfePI/4. ) 

VV2  = VVl*2.*SQRnPI*XX**3)*CEXPH0..1.)*PI/4.)  

VV2=  VV2*4. * ( 0.  , l • ) *CEXP  ( (0. , 1.  ) *P  1/4  • ) *SUR  T(  PI  *XX<‘*3 ) / 3, 

GO  TO  54 

C POWER  SERIES 

53  VV=  Cl.-CEXP(  (0.,l.)*PI/4.r*SQRt(Pt*XX**3)/4. 

l+7.*(0.,i. )*XX**3/60.+7.*SQRT(PI*XX**9)*CEXP< (0.,-l.)*PI/4.)/512.) 

VVl= l.*CEXP( { 0., 1. )*Pl/4.)*SQRT(PI*XX**3)/2.-7.*(0.,l. ) 

1#XX**3/12. -7  #*CEXP ( ( 0. ,-l. )*P 1/4. ) *SQRT ( P I*XX**9 ) / 64. 

W2=l.  + 7.*(0., 1. )*XX**3/12.+7.*CEXP( (0.,-l. )*PI/4. ) *SQRT { P 1*XX*»9 ) 
1/32. 

54  CONTINUE 

Hl=(K**2*Y/(2.*PI))*CEXPU0.,-l.)*K*DIS)/(K*DIS) 

C THE  GREEN'S  FUNCTION 

<HTD»Hl*< ( (0. ,-l. )*COS( T 1**2+ <20. *S I N ( T ) * *2/ 9 •- 29 . *COS IT) **2/ 24 • ) 7 
1 (K*DIS) )*VV) 
v HTS= (0 . ,0.  ) 

RETURN* 

END 
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